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SUMMARY 


This is the final report for NASA Contract NAS8-34427. The report presents the results of 
investigations into the problems of modeling atmospheric backscatter from aerosols, in the 
lowest 20 km of the atmosphere, at CO 2 wavelengths. It contains a summary of the relevant 
aerosol characteristics and their variability, together with a discussion of the measurement 
techniques and the errors involved. The different methods of calculating the aerosol 
backscattering function, both from measured aerosol characteristics, and from optical 
measurements made at other wavelengths, are discussed in detail and limits placed on the ac- 
curacy of these methods. The effects of changing atmospheric humidity and temperature on 
the backscatter are analyzed and related to the actual atmosphere. Finally, the results of 
modeling CO 2 backscatter in the atmosphere are presented and the variation with height and 
geographic location discussed, and limits placed on the magnitude of the backscattering 
function (ico 2 - Detailed conclusions regarding modeling techniques and modeled at- 
mospheric backscatter values are presented in tabular form at the end of the report. The 
most significant of these conclusions are as follows: 

1. The qualitative behavior of fico 2 can be modeled, including the variation with height 
and major geographic features. 

2. Accurate calculation of ftco 2 requires an accurate knowledge of the particle size 
distribution to a particle radius of at least 1 /.<m, for stratospheric aerosols, and at least 5 ^^m 
for tropospheric aerosols. 

3. Particles with radii greather than 0.5 fim are responsible for 50% to 95% of the scat- 
tering from tropospheric aerosols. 

4. Pco 2 decreases from about 10'^ m'* sr' at the Earth’s surface to 10 '* m ' sr* in the 
stratosphere. A minimum value of 10" m * sr * is also observed in the free troposphere over 
the remote ocean. 

5. The variation of pco 2 with wavelength between 9 f^m and 11 (im is not significant, ex- 
cept for water soluble aerosols. 

6. Pco 2 is not strongly dependent upon relative humidity, except for very high values well 
above the point of deliquescence for the water soluble component. 

7. The major information gap is in the composition and size distribution of particles 
with radii greater than 1 ^m, particularly above an altitude of 7 km. 

8. It is recommended that simultaneous measurements be made on aerosols in the free 
troposphere and stratosphere using different particle sensors and CO 2 lidar. 



1. INTRODUCTION 


The recent introduction of coherent doppler lidar systems, for atmospheric wind 
measurements, has led to a need for detailed knowledge of the scattering characteristics of 
stratospheric aerosols, at the CO 2 laser wavelengths used in these experiments. In order to 
meet this need, a program of study has been carried out to assess the problems encountered 
in modeling atmospheric backscatter at a wavelength of 10.6 ^^m (and neighboring 
wavelengths) and to determine, within limits of error, the values of the aerosol backscatter- 
ing function (/3co,) in the atmosphere. 

Factors involved in the calculation of [icoi from measured aerosol characteristics are the 
size distribution and number concentration of the aerosol, its refractive index and the parti- 
cle shape. Consideration must also be given to uncertainties in the measurement of these 
quantities and their natural variation in the atmosphere. These last two factors are par- 
ticularly important and also difficult to assess, as a variety of different measurement tech- 
niques exist, yielding results that often agree qualitatively rather than quantitatively. 
Moreover, such measurements have tended to be made mainly over continental regions, and 
oceanic areas, particularly in the southern hemisphere have been poorly represented. We 
have also examined the possibility of using backscattering measurements made at other 
wavelengths (using visible and near infrared lidar), and extinction measurements, to estimate 
Pcoz- Again, the same factors are found to influence the establishment of any simple rule for 
conversion of backscattering functions and extinction coefficients from one wavelength to 
another. 

Several parallel approaches have been taken to examine these problems. We have used 
log-normal models of variable mode radius and width to simulate the naturally occurring 
aerosol size distributions, using both a flexible range of refractive indices and those cor- 
responding to a set of common aerosol materials. In these and other calculations we have 
assumed spherical particle shape and used standard Mie theory to determine the optical 
properties. We have also carried out similar calculations using a very wide range of 
measured particle size distributions obtained from a comprehensive literature search. 
Calculations of the effects of varying certain atmospheric and optical parameters, e.g., 
humidity and wavelength, have been made using a limited number of measured size distribu- 
tions, representative of certain classes of aerosol measurement. 

The organization of this report is sequential. Section 3 summarizes measured aerosol 
characteristics and the models used to describe them. Variations of concentration with 
altitude and geographical location, as well as the effect of changes in humidity and 
temperature are covered. Section 4 summarizes potentially useful optical measurements 
made on aerosols, using visible and near infrared wavelengths, particularly those made using 
lidar. Sections 5 and 6 discuss the techniques used for the calculation of /Jco^ and the errors 
arising in these techniques. A fairly full and, as far as possible, quantitative discussion is 
given of the influence of the various factors listed above. Section 7 contains the results of 
our calculations on Pco^ the atmosphere, its variation with height and geographical loca- 
tion and its statistical fluctuation. Comparison is made of the rsults of the model calcula- 
tions with available /Icoj lidar data from NASA-Marshall Space Flight Center and from 
NOAA. Section 8 summarizes, in tabular form, our conclusions with regard to the various 
aspects of this work. 



2. DATA BASE USED 


As will be discussed in more detail later in this report, there are three basic methods of 
modeling ftco 2 - The first method is based on our knowledge of the characteristics of aerosol 
particles at all altitudes and geographical locations of interest. The characteristics of aerosol 
particles pertinent to (ico 2 are basically aerosol size distribution, concentration, chemical 
composition, and refractive index. The input data commonly used in this method consists of 
in situ measurments and aerosol particle models. The second and third methods utilize lidar 
data at visible and near infrared wavelengths, and extinction measurements, respectively. 
Although the later two methods are more indirect than the first method, the existence of 
large quantities of data in these two classes makes it worthwhile to investigate the possibility 
of utilizing the lidar data obtained by different research groups scattered around the globe. 
In contrast to the wide distribution of possibly useful lidar data, extinction measurements, 
that serve our purpose, are dominated by the single data set that has been obtained by the 
SAGE and SAM II satellite experiments (McCormick et al., 1979). This data is now in the 
process of being archived at the National Space Sciences Data Center at NASA-Goddard 
Space Flight Center, Greenbelt, Maryland. 

The method of acquisition of data bases, to be used in subsequent studies discussed in this 
report, consisted mainly of literature searches and personal contact. A literature search was 
carried out under the major headings used in this study, using the NASA/RECON system. 
The subject of interest was divided into five major areas which are listed in Table 2. 1 . Details 
of the search procedure as well as the number of references identified are also shown in that 
table. Approximately 1500 sources have been identified in this way. The most important ar- 
ticles have been obtained and analyzed according to their potential value. Besides such 
research reports in open literature, some unpublished data belonging to the Institute for At- 
mospheric Optics and Remote Sensing and its staff scientists have also been examined. 

Our data base has also been augmented by personal contact. This included cor- 
respondence, telephone calls, personal discussions, and a few visits, where possible. This ac- 
tion was confined mainly to the field of lidar studies and details of the results of this in- 
vestigation are shown in detail in Section 4.1. 
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TABLE 2.1 LITERATURE SEARCH ON NASA/RECON SYSTEM 


Subject area 



Details of search 


Number of 

references 

identified 

Set 

No. 

Description 

Aerosol models and 

1 

6546 

6546 

ST/aerosols 


size distributions 

2 

1217 

1217 

BT/particles 



3 

14598 

14598 

R8,R9 ST/aerosols 



4 

2550 

2550 

R10,R13 ST/aerosols 



5 

14 

14 

RT/SAGE satellite 



6 

2876 

2876 

ST/ models 



7 

15502 

15502 

NT/atmospheric models 



8 

599 

599 

NT/environment models 



9 

1917 

1917 

(l+2 + 3 + 4)*(6 + 7 + 8) 



10 

848 

848 

MJ/models 



11 

7846 

7846 

MJ/atmospheric models 



12 

346 

346 

MJ/environment models 



13 

1003 

1003 

(1 +2+3 + 4)*(10+ 11 + 12) 



14 

315 

315 

1*(10+ 11 + 12) 



15 

1396 

1396 

RT/dust 



16 

350 

350 

(l + 2+15)*(10+ll + 12) 



17 

6390 

6390 

R1,R2 ST/size distribution 



18 

1551 

1551 

(1 + 2+ 15)M7 



19 

1380 

1380 

1»17 



20 

2299 

2299 

MJ/R1.R2 ST/size distribution 



21 

695 

695 

(1 + 2+ 15)*20 



22 

670 

670 

21-16 



23 

473 

473 

)22/75-81 

473 

Aerosol composition 

1 

6546 

6546 

ST/aerosols 



2 

1217 

1217 

BT/particles 



15 

1396 

1396 

RT/dust 



28 

744 

744 

R1.R2 ST/composition 



29 

25 

25 

(1 + 2+ 15)*28 



30 

7840 

7840 

ST/chemical composition 



31 

6597 

6597 

ST/chemical analysis 



32 

396 

396 

(l + 2+5)*(30 + 31) 



33 

416 

416 

29+32 

416 

Aerosol refractivity 

1 

6546 

6546 

ST/aerosols 



2 

1217 

1217 

BT/particles 



10 

848 

848 

MJ/models 



11 

7846 

7846 

MJ/atmospheric models 



12 

346 

346 

MJ/environment models 



15 

1396 

1396 

RT/dust 



16 

350 

350 

(l + 2+15)*(10+ll + 12) 



20 

2299 

2299 

MJ/R1.R2 ST/size distribution 



21 

695 

695 

(1+2+ 15)*20 



22 

670 

670 

21-16 



23 

473 

473 

)22/75-81 



24 

4771 

4771 

US/refractivity 



25 

233 

233 

(1 +2 + 15)‘24 



26 

220 

220 

25-16 



27 

188 

188 

25-23 

188 

Lidar measurements 

1 

2276 

2276 

us/optical radar 


of atmospheric 

2 

6546 

6546 

ST/aerosols 


aerosols 

3 

1217 

1217 

BT/particles 



4 

1396 

1396 

RT/dust 



5 

440 

440 

l*(2 + 3 + 4) 

440 

Extinction radiometry 

1 

409 

409 

ST/extinction 


and optical depth 

2 

4500 

4500 

ST/atmospheric attenuat 



3 

2988 

2988 

ST/atmospheric boundary 



4 

1763 

1763 

MJ/atmospheric optics 



5 

3844 

3844 

ST/ radiometers 



6 

24 

24 

4*5 



7 

88 

88 

l*(2+3 + 4) 



8 

42 

42 

UTP/optical *3 depth 



9 

644 

810 

AX/optical *3 depth 



10 

80 

115 

(2 + 3 + 4)*(8 + 9) 



11 

184 

219 

6+7+10 

184 


Terms used in the search are combined together using the normal rules of Boolean algebra. Special symbols usesd are as follows: 


AX 

Key word occurring in abstract 

RT 

Related Term 

BT 

Broader Term 

ST 

Subject Term 

MJ 

Major Term 

US 

Use 

NT 

Narrower Term 

UTP 

Key word occurring in title 
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3. AEROSOL PROPERTIES— MODELS AND MEASUREMENTS 


3.1 Survey 

Aerosol properties such as composition, refractive index and size distribution are a func- 
tion, not only of geographic location and height above the earth’s surface, but also of season 
and the history of the air-mass in which the aerosol occurs. A simplified view of its depen- 
dance upon the first two of these parameters is shown in Fig. 3.1. In this figure the first 20 
km of the earth’s atmosphere has been divided into four main regions within which the 
variation of aerosol properties is less than it is for the whole globe. 

Two classes of boundary layer aerosols are recognized, a maritime aerosol with a moist 
sea-salt aerosol as its main component and a continental aerosol consisting of a mixture of 
an insoluble component (mainly silicate dust and soil) and a soluble component (mainly am- 
monium sulphate and sulphuric acid)(WMO report No. WGP-12).* Exchange of material 
between these two regions will occur near the continental margins. Above the boundary 
layer, there exists a free tropospheric aerosol whose constituents are derived from, and are 
the same as, the continental boundary layer aerosol. This is approximately true whether we 
consider the free troposphere over the land or the ocean although, as we shall see later, the 
size distributions are different. 

Above the troposphere there exists a .background stratospheric aerosol consisting of a 
mixture of approximately 75% sulphuric acid and 25% water. This is distributed globally 
with a concentration that varies with height and latitude as shown in Fig. 3.2 (Rosen et al., 
1975). 

Also shown in Fig. 3.1 is the approximate height variation of the aerosol number concen- 
tration. This quantity is highly variable but, in general, exhibits two maxima, one near the 
surface of the earth and the other in the stratosphere. The concentration near the earth’s sur- 
face depends on the amount of mechanical disturbance; that in the stratosphere depends on 
the recent history of volcanic eruptions which inject material into that region. The concen- 
tration in the free troposphere is also highly variable, particularly over the oceans, where it 
depends upon transport of material from over the continents. 

Fig. 3.3 shows typical examples of the aerosol size distribution in different parts of the at- 
mosphere. These have been plotted in two ways — as the number of particles per log radius 
interval (dN/d log (r)) and as the particle volume per log radius interval (dV/d log (r)). Both 
representations occur in the literature, the former being the more common and the latter 
more useful for conversion to mass aerosol loading. A significant feature of the boundary 
layer and free tropospheric distributions is their bimodal character (Patterson and Gillette, 
1977, Patterson et al., 1980). Of most importance is the coarse particle mode. This consists 
of approximately 80% insoluble and 20% soluble material, has a volume mode maximum 
radius between 1 and 4 /im (in the free troposphere) and a minimum between it and the ac- 

*Note: A third class of Urban-Industrial Aerosols is also recognized, which will not be con- 
sidered in this report. 
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cumulation mode aerosol at about 0.2/im (Nilsson, 1979; Winkler, 1974). The stratospheric 
aerosol is usually considered to be uni-modal in character although recent measurements by 
Hofmann and Rosen (1981) have cast some doubt upon this. 

Measurements on the size distribution of the boundary layer aerosol are well distributed 
over the globe; those in the free troposphere and stratosphere are much less so. Fig. 3.4 
shows the locations of those free tropospheric and stratospheric aerosol measurements used 
in this report. Three continents only are represented. North America, Australia and Europe, 
and a limited amount of data is available over the Atlantic and Pacific oceans. 


3.2 Size Distributions and Concentrations 
3.2.1 In situ measurement techniques 

The existence of a wide range of aerosol size distributions, concentrations and refractive 
indices presents formidable problems in the measurement of aerosol size spectra. The size 
range of aerosol particles in the atmosphere that may contribute to the backscatter at CO^ 
wavelengths extends from about 0.01 /im to 30 jxm and the overall concentration may vary 
from less than 1 cm ’ in the stratosphere to 10’ cm ’ in a polluted area. Aerosol refractive in- 
dices vary from 1.18 - 0.0677i (for water at A = 10.6 f^m) to 2.22 - 0.73i (for soot at A = 10.6 
/im). The problem is further complicated by the fact that most measurement devices perturb 
the natural state of the atmospheric aerosols in advance of measuring the number concentra- 
tion. Some particles may be lost on their passage through the narrow tube leading to the 
counter and the size and composition may be modified deliberately or inadvertently inside 
the instrument. Though theoretically the actual aerosol size distribution can be deduced 
from the measured result, in practice the conversion process depends on so many factors 
that different aerosol size distributions may be obtained from the same set of measurements. 
Consequently, the reported aerosol size distribution obtained from in situ measurements 
must be used with caution and the possible errors involved in the measurement should be 
considered carefully. 

According to Lundgren (1976), aerosol measurement instruments can be divided into four 
categories; inertial impactor, light-scattering particle counters, electric mobile analyzers and 
condensation nuclei counters. A critical discussion of instruments in these four categories is 
given in the following paragraphs. 

A. Inertial impactor 

Inertial classification devices separate aerosol particles in an aerodynamic size dependent 
manner for subsequent analysis by other techniques. The measured aerodynamic size range 
is from about 0.25 /<m to 2.5 /^m in radius. Some possible sources of errors commonly 
associated with all inertial impactors are: 

(a) Aerosols collected by this method may undergo evaporation, growth and other 
physical or chemical processes between the time and location of collection and the time and 
location of analysis; i.e., this method is not a truly in situ method. 

(b) Aerosol particles may be lost in the walls of inlet system; i.e., isokinetic sampling can 
hardly be achieved. 

(c) Since the aerosol size distribution is not obtained directly from the impactor device, 
different size distributions may be deduced from the same set of samples using different 
analyzing techniques. 
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B. Light scattering particle counters 

The basic parameter measured in this class of devices is the amount of radiation scattered 
by the aerosol particle inside the counter. The number and height of the pulses generated by 
the scattered radiation is electrically converted to a particle number and particle size; most 
optical counters operate best for particles with radii between 0.25 p^m and 2.5 ^m. Some 
possible sources of error commonly associated with all optical particle counters are: 

(a) Systematic errors in the electric and optical system such as the instability of the light 
source, lens system distortion and contamination of the optical system. 

(b) Uncertainty of the properties of aerosol, including its shape and the wavelength 
dependence of its refractive index. 

(c) Multi-valued light scattering versus particle size response and multi-particle scatter- 
ing. 

C. Electric mobile analyzer 

This device measures the electric mobility of an unknown particle in a known electric 
field. Since the electric mobility of a gas borne particle in an electric field is a monotonfc 
function of particle size, the size of aerosol particles can be determined from the measured 
electric mobility. There are three basic parts of an electric mobile analyzer: a diffusion 
charge section which places a charge on each aerosol particle; a mobility analyzer which col- 
lects particles with a mobility greater than a pre-set value; and an electrometer which 
measures the current flow by the remaining charged particles. 

This method is applicable to the measurement of aerosol particles with radii less than 
about 0.25 /um. In order to measure larger sized aerosols, it is necessary, either to increase 
the field intensity in the condenser, or to reduce the rate of flow of the aerosol. Very intense 
fields may result in leakages, local discharge and, at low flow rates, other complications such 
as the generation of convection currents. 

Other possible sources of error are the non-uniform charging of the aerosol particles and 
non-laminar flow through the condenser. 

D. Condensation nuclei counter 

The condensation nuclei counter (C.N.C.) is a device that counts the number of active 
nuclei in a supersaturated chamber. The counting can be done manually, or automatically, 
by measuring the extinction of light caused by the presence of liquid droplets containing 
condensation nuclei. 

Two types of errors may occur. One is the instrumental error and another is caused by the 
nature of the particle being measured. The traditional technique for detecting condensation 
nuclei is to make them grow to sizes where they are efficient scatterers of light, by super- 
saturating the sampling air. Since the number of active scattering nuclei increases 
dramatically in response to a slight increase of supersaturation, the supersaturation in the 
chamber must be accurately determined. Furthermore, the growth of condensation nuclei 
depends on their physical characteristics, including composition and shape. If the concentra- 
tion is deduced by measuring extinction, one additional parameter, namely the index of 
refraction is involved. If the characteristics of the in situ aerosol particles are quite different 
from those used for a calibration, large errors may be produced. 

A very brief summary of the common in situ measurement techniques is given in Table 
3.1. The list of names of users shown in the Table is of course incomplete. Most of the 
names listed are here because they have measured aerosol characteristics (particularly the 
aerosol size distribution), which are specifically discussed in this report. 


3-3 



3.2.2 Measured size distributions 


Aerosol size distributions in the free troposphere show considerable differences when 
measurements made by independent workers, using different techniques, are compared. Fig. 
3.5 shows a set of volume size distributions obtained by Blifford and Ringer (1969) at dif- 
ferent altitudes over the continental U.S.A. using a single stage impactor. All the distribu- 
tions show a mode maximum at a radius of about 2 pm. There is a decrease of concentration 
with altitude particularly for the smaller particles. Fig. 3.6 shows size distributions measured 
by Cress (1980) over Europe using a Royco optical counter. Although not exactly the same, 
the general characteristics of the curves are similar, the mode maxima for this set occurring 
around 3 pm. Fig. 3.7 shows a set of distributions obtained by Patterson et al. (1980) over 
the continental U.S.A. and the Pacific Ocean using a PMS optical counter. The boundary 
layer distributions have a shape similar to those obtained by Blifford and Ringer and by 
Cress but the free tropospheric distributions show significant differences. From the point of 
view of modelling the most important difference is that of the mode maximum for the 
large particle distribution, which is less than 1 pm for this data. The effects of these dif- 
ferences on the modelled values for /Jcoj will be discussed in a later section. 

In Section 3.1 we noted that the free tropospheric aerosol over the ocean was largely deriv- 
ed from that over the continent. This has important implications with regard to the size 
distribution due to the gradual sedimentation of the larger particles during transport over 
the ocean. This is illustrated in Fig. 3.8 which shows data obtained by three independent sets 
of research workers over the land and over the ocean. That of Savoie (1978) is related to the 
movement of Saharan dust over the Atlantic Ocean and the change in the mode radius bet- 
ween the size distribution over Sal Island (off the coast of West Africa) and Miami is very 
clear. The other two sets do not relate directly to a specific transport process but illustrate 
the general differences between the aerosol over the Pacific Ocean and the continental 
U.S.A. In both cases the free troposphere over the ocean contains far fewer large particles. 
It is also clear from recent work (Shaw, 1980; Duce et al., 1980) that desert dust is carried 
over the Pacific Ocean from Eastern Asia in a somewhat similar manner to that in which 
Saharan dust is carried over the Atlantic. At this time measurements on these dust streams 
have yet to be made in the free troposphere. 

An alternative way to present the differences between the free troposphere over the ocean 
and the land is to show the height variation of the larger particles for both cases. This is 
shown in Figs. 3.9, 3.10 and 3.11 for data over Europe, U.S.A. and the Pacific Ocean, 
respectively. The rapid decrease in concentration above the marine boundary layer is most 
noticeable. This is in contrast to the continental cases where convection causes vertical mix- 
ing throughout the troposphere. 

Patterson et al. (1980) have reported a strong latitudinal variation in the aerosol mass 
loading in the free troposphere over the Pacific Ocean. The variation, which is shown in Fig. 
3.12, reflects the fact that the majority of the continental land mass and hence the most 
significant movement of continental dust over the oceans is in the Northern Hemisphere. 
(The scale on the right hand side of this figure, which shows the equivalent value of /Jcoj will 
be discussed in Section 7.) 

So far in this section we have been primarily concerned with the variation with height and 
geographical location of the mean value of the aerosol concentration. It is of importance to 
inquire about the variation that may be observed if these two parameters are held fixed. The 
same data sets, as previously used from Cress and Patterson et al., have been analyzed in 
order to obtain information on this variation. Fig. 3.13 shows the probability distributions 
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for the particle count obtained by Cress at three different altitude ranges for a nine month 
period over Europe. No division in particle size range has been attempted although the data 
indicate that there was slightly more variation for particles with radii greater than one 
micron than for particles smaller than this. Fig. 3.14 shows probability distributions 
obtained by Patterson et al. on flights at a height of 5-6 km over the western U.S.A. and the 
central Pacific. Although the data are somewhat limited, the range of statistical variation is 
similar to that obtained by Cress, namely that the 10% and 90% probability levels for the 
particle concentration lie approximately 1 Vi decades apart. 

Background stratospheric aerosol size distributions are traditionally uni-modal (Russell et 
al., 1981) and approximately log-normal in shape with a maximum (in dN/d log r) near a 
radius of 0. 1 The concentration shows a latitudinal variation as already described in Fig. 
3.2 and small seasonal fluctuations. The major factor controlling the concentration of 
stratospheric aerosols is, however, volcanic activity, which injects SO 2 (and ash) into the 
stratosphere. This is subsequently oxidized and combines with water to form H 2 SO 4 which 
then associates to form aerosols. The lifetime of these aerosols is to be measured in months 
or years and the injection from a single volcano may eventually spread to cover the globe. 
The ash which is injected with the SO 2 consists of larger particles which sedimate out much 
more rapidly. Although stratospheric particle concentrations have been measured directly 
for the past twenty years, one of the best long-term records is that obtained using lidar 
techniques. This will be described later in Section 4.1. Recently, as noted in Section 3.1, 
Hofmann and Rosen (1981) have made measurements, using a specially designed large parti- 
cle counter, of the concentrations of particles with radii between 1 and 2 /.tm. An example of 
these measurements is shown in Fig. 3.15, where it can be seen that, compared to the log- 
normal distribution, there is a significant excess of particles with radii within these limits. 
Unfortunately, no other workers have yet made any similar measurements and the nature 
and behavior of these particles still has to be established. 

In subsequent sections of this report we shall make use of both number and mass concen- 
trations of aerosols when modelling (icoi- order to assist in the appreciation of those sec- 
tions, Table 3.2 shows a brief list of typical values of these quantities for different altitudes 
and parts of the globe. A summary of the main features of aerosol size distributions and 
concentrations, most of which have been covered in this chapter, is given in Table 3.3. 


3.3 Aerosol Models 

In the past two decades, numerous in situ measurements of aerosol particles in the at- 
mosphere have been reported. These experimental data show that, in general, the physical, 
optical and chemical properties of aerosol particles change with altitude, location and season 
of the year. In order to predict the behavior of aerosol particles in different regions of the at- 
mosphere and under different ambient conditions, several investigators have summarized 
the reported aerosol measurements and derived analytical or numerical representations of 
aerosol properties. A survey of these global aerosol models shows that their number is sur- 
prisingly small (see Table 3.4). Furthermore, most of these investigators have concentrated 
their effort on the lower stratosphere and boundary layer; the free troposphere is the region 
with the least amount of information available. This is probably due to the fact that aerosol 
particles in the free troposphere are more difficult to measure than those in the boundary 
layer, and the variation of aerosol properties in the free troposphere is much larger than in 
the lower stratosphere. We have also found that most of the reported aerosol models consist 
of information on the size distribution and refractive index but not on the number concen- 
tration. On the other hand, even though aerosol size distributions in the stratosphere are 
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meticulously simulated, the size ranges of interest are usually for the smaller particles, the 
occasionally reported existence of larger particles (> 1 f^m) in the stratosphere has never been 
included in the stratospheric models. In summary, we found that the properties of aerosol 
particles in the free troposphere, the number concentration of aerosol particles, and the size 
distribution of larger particles in the stratosphere are three areas of information not ade- 
quately simulated in the reported aerosol models. 

Some of the global aerosol models important for /Jcoj studies are summarized in Table 3.4. 
As indicated earlier the atmosphere of interest can be divided into three regions; lower 
stratosphere (from tropopause to 30 km), the free troposphere (from '^2 km to tropopause), 
and boundary layer (from ground to ~2 km). Two types of aerosols in the stratosphere are 
modeled— background and post-volcanic. Aerosols in the boundary layer can be summa- 
rized in three categories: maritime, clean air continental, and polluted continental. 

It should be noted that aerosol models listed in Table 3.4 are by no means inclusive. Some 
models, such as the Junge power law and the Deirmendjian modified gamma are not listed 
separately because these have been utilized and incorporated into the aerosol models 
developed by later investigators. It may also be noted that many models are based on ex- 
perimental data obtained over limited regions of the earth’s surface, e.g., the stratospheric 
model of Toon and Pollack is based only on data from the USA and Australia. Care must be 
taken in applying such models to conditions outside those from which they have been 
developed. 

3.4 Composition and Refractive Index 

As mentioned in Section 3.1 of this report, the composition of aerosol particles depends 
on several parameters, including the geographical location, height, season and history of the 
air mass containing the aerosol particles. Roughly speaking, the variation of aerosol com- 
position decreases with height. In the boundary layer, the continental aerosols consist of a 
mixture of soluble components (mainly ammonium sulfate and sulfuric acid) and an insolu- 
ble component (mostly silicate dust and soil); maritime aerosols are mainly sea-salt particles 
but a continental component may also be present. Above the boundary layer, but below the 
tropopause (free troposphere), the composition of the aerosol particles is similar to that in 
the continental boundary layer. Above the tropopause, the aerosol particles are mainly 
sulfuric acid solutions (approximately 75 %H 2 S 04 and 25'%H20). 

One of the important input parameters in calculating ftcoi from the Mie code is the com- 
plex refractive index of the aerosol at the wavelength under consideration. The wavelength 
dependence of the complex refractive indices of some common aerosol materials are listed in 
Tables 3.5 and 3.6. In Table 3.5 most of the refractive indices are for visible or near infrared 
wavelengths. In addition, the refractive indices at 5.0 ^m and 10.6 ^m are also listed for the 
sake of comparison. In Table 3.6 the wavelength of radiation is limited to CO 2 laser 
wavelengths, ranging from 9.1 fim to ll.l f.»m. It should be kept in mind that for some 
substances there are slight discrepancies between refractive indices given by different in- 
vestigators. When this is the case, we have chosen one representative set of values to be listed 
in the tables. Furthermore, some of the limited values are obtained by interpolation and 
some are obtained by very rough estimation from graphs published by the authors. It should 
also be noted that the refractive indices listed (except water) are those of dry aerosols at 
room temperatures. A change of the ambient parameters such as relative humidity and 
temperature wiil certainly affect the refractive index. The influence of these on the refractive 
index will be examined in the next subsection. 


3-6 



From the listed refractive indices shown in Tables 3.5 and 3.6, the following conclusions 
can be drawn; 

(a) In general, the refractive indices are quite sensitive to CO 2 wavelengths but not sen- 
sitive to visible or near infrared wavelengths. 

(b) In general, the imaginary parts of the refractive indices at visible wavelengths are so 
small that they can be ignored in Mie calculations; however, the imaginary parts of the 
refractive indices at COj wavelengths, for most materials except NaCl, are quite large and 
the absorption by aerosols at these wavelengths should not be ignored in any Mie calcula- 
tion. 

(c) At visible wavelengths the variation of refractive index with wavelength is monotonic, 
but at CO 2 wavelengths the variation is not monotonic for some materials. Consequently, in- 
terpolation of the refractive indices at CO 2 wavelengths must be carried out with caution. 

(d) At CO 2 wavelengths, the refractive indices of common aerosol materials are, in 
general, larger than those at visible wavelengths. Furthermore, for some aerosol materials, 
the refractive index (at certain wavelengths) may be much larger than at neighboring 
wavelengths [e.g., for (NH 4 ) 2 SOi at wavelength A = 9.5 ptm]. The backscatter at such a 
wavelength may be quite different from that at other wavelengths. 


3.5 Humidity and Temperature Effects 


As the relative humidity increases, water vapor may be absorbed by aerosol particles 
suspended in the atmosphere. Laboratory experiments and theoretical calculations have 
shown that the properties of aerosol particles do not change until the ambiednt humidity has 
reached a threshold value which depends on the deliquescent property of the aerosol com- 
pound under consideration. The condensed water will change the optical properties in two 
ways: the first is to increase the size of the aerosol particles and the second is to change their 
effective refractive index. 

The variation of the size of aerosol particles with ambient humidity can be modeled in 
several ways. The easiest way is to utilize the published growth factors which take into ac- 
count the change of aerosol size with relative humidity. Unfortunately, most of the pub- 
lished data on growth factors are limited to certain types of aerosols in the boundary layer. 
For pure hygroscopic compounds in the atmosphere, the growth factor can be calculated as 
follows. 


Since the mass of solute in an aerosol droplet before and after changes of ambient humidi- 
ty should be the same, we have. 

vex = v'e'x' (1) 

where v is the volume of one aerosol particle, q its density and x the weight percentage of 
solute in the solution. The symbols without and with prime refer to the initial, and final con- 
dition, respectively. 


From Eq. (1) we have 

V' _ QX 
V ~ 


( 2 ) 


and for the growth factor, defined as the ratio of radii after and before the change of am- 
bient conditions. 



r 



( 3 ) 
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The new density q ' is a function of x' and the new composition x' can be estimated by con- 
sidering the fact that the ambient water vapor pressure should be in equilibrium with the 
water vapor pressure above the surface of the droplet solution with composition x'. 

As mentioned earlier, the change of ambient humidity may change the aerosol composi- 
tion, resulting in a change in its effective refractive index. The variation of refractive index 
with relative humidity is given by the following equation (Hanel, 1976). 

n = n„ -t- (n„ - nj -t- A + — • (4) 

V 6w 

where n is the real or imaginary part of the complex index of refraction of the aerosols at the 
relative humidity under consideration, n, q , and m with subscripts refer to the real or imagi- 
nary part of the complex index of refraction, density and mass, respectively, subscript w 
refers to water and subscript o refers to the pure substance at 097o relative humidity. The 
values of m,, and m„ are functions of the ambient relative humidity and they can be deter- 
mined once the value of x' has been found with the method described earlier. 

So far we have discussed only the effect of humidity on aerosol size and refractive index. 
If we assume the water content of the atmosphere remains constant, the relative humidity 
will change following a change in ambient temperature. Consequently, due to the change in 
relative humidity, aerosol size and refractive index should vary with the ambient 
temperature. Additionally, there may be some direct changes in the aerosol density and 
refractive index with temperature, though such changes are generally small. The change of 
densities of pure substances are given by standard references such as the International 
Critical Tables and the effect of temperature on the index of refraction can be derived from 
the well-known Lorentz-Lorenz formula (Longhurst, 1964): 


n^ - 1 


= constant 


( n ^ + 2)q 

The value of n at temperature T is then given by 


SiT) + 1 

- [ 1 - CeCnJ 


‘/2 


(5) 


( 6 ) 


Where the constant C can be calculated from the known index of refraction n' at 
temperature T' 


C = 


n'HT') - 1 
[n'HT') -t- 2] e(T') 


(7) 
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TABLE 3.1: IN SITU 

INSTRUMENTS USED TO I4EASURE AEROSOLS IN THE 

ATMOSPHERE 


TYPE 

NAME 

USERS 


INERTIAL IMPACTOR 

FILTER COLLECTION 
PARTICLE IMPACTOR ON SLIDES 

PARTICLE IMPACTOR ON SCREENS 
WIRE COLLECTOR 

QUARTZ CRYSTAL MICROBALANCE (QCM) 

SAVOIE (1978), PROSPERO (1972) 
BLIFFORD& RINGER (1969) , DELUISI 
et al (1976), HOBBS et al (1976) 
GRAS AND AYERS (1979) , BIGG (1980) 
FARLOW et al (1981) 

WOODS (1979) 

LIGHT SCATTERING 
PARTICLE COUNTER 

1 

ROYCO OPTICAL COUNTER 

SINGLE PARTICLE OPTICAL COUNTER 

DUSTS ONDE 

KNOLLENBERG PROBE 

INTEGRATED NEPHELOMETER 

POLAR NEPHELOMETER 

CRESS (1980) 

PATTERSON (1980) 

HOFMANN and ROSEN (1981) 
HOBBS et al (1976) 

CRESS (1980) 

GRAMS & PATTERSON (1980) 


ELECTRIC MOBILE ANALYZER 

ELECTRIC AEROSOL ANALYZER (EAA) 

HOBBS et al (1976) 


CONDENSATION NUCLEI COUNTER 

GARDNER SMALL PARTICLE COUNTER 
CCN COUNTER 
ICE NUCLEI COUNTER 
STRATOSPHERIC AITKEN NUCLEI 
DETECTION SYSTEM (SANDS) 

HOBBS et al (1976) 
PUESCHEL (1980) 
HOBBS et al (1976) 

HABERL (1975) 




TABLE 3.2: TYPICAL AEROSOL CONCENTRATIONS 


NUMBER CONCENTRATIONS 


RADIUS, r N(>r), cm 


TYPE 

AUTHOR 




STRATOSPHERIC, 

ROSEN & HOFMANN 

.25 


.25 

PRE-VOLCANIC 

(1975) 




STRATOSPHERIC, 

GRAS (1976) 

.22 


o 

I — 1 

POST-VOLCANIC 





CONTINENTAL, 5 KM 

BLIFFORD (1970) 

.13 

Hm 

'V^l.O 

MARINE , 5 KM 

BLIFFORD (1970) 

.13 

nm 

'^^ O.l 

CONTINENTAL, 5 KM 

PATTERSON et al (1980) 

.5 

nm 

''^ 0.3 

CONTINENTAL, 0 6 KM 

CRESS (1980) 

.4 

fim 

0.1 — >1 


MASS CONCENTRATION 



-3 

TYPE 

AUTHOR 



M,/ig m 

MARINE , 5 KM 

PATTERSON et al (1980) 



. 02 — 1.0 

CONTINENTAL, 5 KM 

PATTERSON et al (1980) 



'^ 1.0 

ATLANTIC DUST, 3 KM 

PROSPERO & CARLSON (19721 



50 

ATLANTIC DUST, 0 KM 

PROSPERO & CARLSON (1972) 



1 ^ z u 

PACIFIC DUST, 0 KM 

DUCE et al (1980) 



1 —2 

PACIFIC DUST, 2-3 KM 

SHAW (1980) 



0,50 

SOUTH POLE, 0 KM 

BIGG (1979) 



> 0.14 
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TABLE 3.3: SUMMARY OF THE MAIN PROPERTIES OF AEROSOLS AND 

THEIR DISTRIBUTION 


1. GEOGRAPHICAL DISTRIBUTIONS ARE HARD TO DETECT BUT 
SIGNIFICANT DIFFERENCES EXIST BETWEEN THE FREE TROPO- 
SPHERIC AEROSOL OVER THE LAND AND THE OCEAN. 

2. THE FREE TROPOSPHERIC AEROSOL OVER THE OCEAN IS HIGHLY 
VARIABLE IN CERTAIN REGIONS DUE TO TRANSPORT OF CONTINENTAL 
DUST. 

3. THE PRINCIPAL VARIATION IN AEROSOL CONCENTRATION IS THAT 
WITH ALTITUDE. IT IS DIFFERENT OVER LAND AND OCEAN. THE 
SMALL PARTICLE CONCENTRATION FALLS MORE RAPIDLY THAN THE 
LARGE PARTICLE CONCENTRATION 

THE CONCENTRATION OF STRATOSPHERIC AEROSOL IS RELATED TO 
VOLCANIC ACTIVITY^ REACHING A PEAK SHORTLY AFTER A MAJOR 
ERUPTION AND FALLING QUAS I “EXPONENT I ALLY WITH A TIME 
CONSTANT OF 6“12 MONTHS. 

5. THE UPPER TROPOSPHERE IS NOT WELL DOCUMENTED BUT MATERIAL 
IS PROBABLY EXCHANGED THROUGH THIS REGION BETWEEN THE 
STRATOSPHERE AND LOWER TROPOSPHERE. 

6. SIGNIFICANT VERTICAL STRATIFICATION AND HORIZONTAL VARIATION 
EXISTS. 

7. DIFFERENT MEASUREMENT TECHNIQUES YIELD SIGNIFICANTLY DIFFERENT 
VALUES FOR DN/d(L0G r) 
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TABLE 3,4: Summary of Global Aerosol Models Important for Studies 


Report of the 
meeting of JSC 
experts (Geneva, 
27-31 Oct 1980) 
Chairman: H.J. 
Bolle 


Hanel & Bullrich 


Kuznetsov and 
Izhovkina 


4 I McClatchy et al 


Aerosols and Climate 1980k.tratof.pher 


free 

troposphere 


Physical Properties 
of Four Types of 
Tropospheric Aerosols 


planetary 

boundary 

layer 


lower 

troposphere 


Two Models of the 
Atmospheric Aerosol 


whole 


Optical Properties 
of the Atmosphere 
(third edition) 


1972 whole 

atmosphere 
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(Continued) 


No. 

5 


6 


7 


8 


Author 

Pinnick et al. 


Russell et al. 


Shettle & Fenn 


Toon and Pollack 


Title 

Year 

Applicable 

Regions 

Height Range 
& resolution 

Stratospheric Aerosol 
Measurements IV. 
Optical Model 
Calculation 

1976 

stratosphere 

18 20 km 

Satellite and 
Correlative Measure- 
ments of the Strato- 
spheric Aerosol. I. 

1981 

stratosphere 

T+2 30 km 

An Optical Model for 
Data Conversions 


tropopause 

troposphere 

0.7T -*• T+2 

0.3T -*• 0.7T 

T = tropopause 
height 

Models of the Aerosol: 

1979 

troposphere 

above boundary 

of the Lower Atmo- 


lower 

layer 

boundary layer 

sphere and the Effect: 
of Humidity Varia- 
tions on Their Opti- 
cal Properties, 
AFGL-TR-79-0214 


altitudes 

A Global Average 
Model of Atmospheric 
Aerosols for Radio- 
active Transfer 
Calculations 

1976 

stratosphere 

12 -*■ 45 km 
[resolution 3 km) 

0 ■+ 12 km 
(resolution 3 km) 


Environmental 

Conditions 


over Laramie 
1971 - 1974 


a) background 

b) post volcanip 
eruption 
global avg. 

global average 

global average 


a) rural 

b) urban 

c) maritime 


a) background 

b) post-volcani)c 
eruption 

global average 


Comments 


analytic aerosol 
size distribution 
and refractive 
index based on 
own measurements 

summarizes models 
used previously 


size distribution 
is the sum of two 
log-normal distri- 
butions 


size distributions 
are Zold and power 
law 
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TABLE 3.5: 


REFRACTIVE INDICES OF CO»«»l AEROSOL MATERIALS AT COj WAVELENGTHS 






COMPLEX INDEX 

OF REFRACTION 









wavelength in microns 




AUTHORS 

MATERIAL 

9.1 

9-2 

9.5 

9.8 

10.0 

10.59 

11.0 

U.l 


Water 

1.259- 0. U05 

1.255- 0- 0415 

1. 243 - 0.0444 

1.229- 0.0479 

1.218- 0.0508 

1.179-0.0723 

1.153 - 0.0968 

1.138-0. 106 

Hale & Querry 
(1973) 

Water 

soluble 

2. 330- 0.395 

2.200- 0.420 

1.950 - 0.160 

1.870 - 0,095 

1.820- 0.09 

1.760- 0.07 

1.720 - 0.05 

1.710-0.049 

Shettla £ r«nn 
(1979) 

Dust- 

iike 

1.710- 0.145 

1,720- 0.150 

1.730- 0.162 

1.740 - 0. 162 

1.750* 0.162 

1.620-0.120 

1.620-0.105 

1.614 - 0.104 

Shettle & Perm 

(19793 

Soot 

2.175- 0.^00 

2.180 - 0.700 

2.190 - 0.710 

2.200 - 0.715 

2.210- 0.720 

2.220- 0.730 

2.230 - 0.730 

2-232 - 0.732 

Shettle fi Penn 
(1979) 

Sea Sait 

1.630 - 0.027 

1.610 - 0.026 

1. 580 - 0.018 

1.560- 0.016 

1.540- 0.015 

1.500- 0.0X4 

1.480 - 0.014 

1.480- 0.014 

Shettle & Penn 
(1979) 

clay 

- O. 1 

- 0. 35 

- o.r.5 

- 0.6 

- 0.5 

- 0.2 

« 0. 13 

- 0.16 

Patterson 

(1981) 

7S\ 

1.651 - 0.‘;54 

1.615 - 0,546 

1.631 - 0.699 

1.944 - 0.538 

1.915-0.396 

1.737 - 0.273 

1.676- 0.410 

1.739 - 0.463 

Palmer & Williams 
(1975) 


1.172- 0.0093 

1.U4- 0.010 

1.06 -0.013 

0.976 - 1.79B 

Q.92 -0.018 

0.55 -0.06 

0.30 - 0.089 

0.275- 0.179 

Toon et ai. 
(1976) 

NaC« 

l.J99-C.l::10‘® 

1.498-0. Ixl0‘® 

1.495-0. IxIO"® 

l. 492-0. Ixlo’® 

1.49 - O.lxio'^ 

1.49 - O.lxio'® 

1,49 - 1.9xl0”’ 

1.489-1.9x10"’ 

Toon et al. 

U976) 

(NH4)2S0^ 

i.-;6 - 1.9 

1.95 - 2.1 

2.68 -0.53 

2.376- 0. 24 

2.19 -0.13 

1.972 - 0.057 

1.90 - 0.04 3 

1.88 -0-04 

Toon et dl. 
(1976) 

Sahara 

dust 

1. JO - 0.6 

2.0 -0.65 

2.9 - 0.75 

2.8 -0,85 

2.7 -0.9 

1.74 -0.60 

1.85 -o.je 

i.ae -0.36 

Volz (1973) 

Volcanic 

Jui^t 

1.00 - 0.: 

2.Q -0.2S 

3.2 -0.35 

3.1 -0.35 

3.0 -0.4 

1.95 -0.40 

1.90 - 0. 38 

1.88 -0,36 

Volz (1973) 
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TABLF: 3.6: 


REFRACTIVE INDICES OF COMMON AEROSOL MATERIALS 


AEROSOL 



COMPLEX 

INDEX OF REFRACTION 




MODELS 

0.488 

0.55 

wavelength in microns 
0.6328 1.06 

1.66 

5.0 

10.59 

AUTHORS 

Water 

1.336- Ixio"® 

1.333- 1.96X lo"® 

1.332 - 1.5x 10*® 

1.326- 5x10*® 

1.316- 9.43 x 10*® 

1.325- 0.0124 

1.179 - 0.0677 

Hale 6 
Querry 1973 

Maritime 

1.418 - 0.002 

1.418 - 0.002 

1.415 - 0.002 

1.405 - 0.004 

1.376 - 0.004 

1.372 - 0.010 

1.380 - 0.057 

Selby 1976 

Maritime/ 

Rural 

1.475 - 0.005 

1.474 - 0.004 

1.473 - 0.004 

1.463 - 0.009 

1.408 - 0.010 

1.381 - 0.012 

1.550 - 0.071 

Selby 1976 

Rural 

1.530 - 0.008 

1.530 - 0.006 

1.530 - 0.006 

1.520 - 0.014 

1.440 - 0.016 

1.390 - 0.013 

1.720 - 0.085 

Selby 1976 

Rural/ 

Urban 

1.569 - 0.086 

1.569 - 0.082 

1.569 - 0.080 

1.560 - 0.089 

1.500 - 0.096 

1.492 - 0.116 

1.810 - 0.198 

Selby 1976 

Urban 

1.607 - 0.163 

1.607 - 0.158 

1.607 - 0.154 

1.600 - 0.163 

1.559 - 0.175 

1.593 - 0.218 

1.895 - 0.310 

Selby 1976 

Water 

Soluble 

1.530 - 0.005 

1.530 - 0.006 

1.530 - 0.006 

1.520 - 0.017 

1.487 - 2 X 10*’ 

1.450 - 0.012 

1.760 - 0.07 

Shettle & 
Fenn 1979 

Dust- 

like 

1.530 - 0-008 

1.530 - 0.008 

1.530 - 0.008 

1.520 - 0.008 

1.367- Bx 10*’ 

1.250 - 0.016 

1.620 - 0.120 

Shettle & 
Fenn 1979 

Soot 

1.750 - 0.450 

1.750 - 0.440 

1.750 - 0.430 

1.750 - 0.440 

1.78 - 0.469 

1.970 - 0.600 

2.220 - 0.730 

Shettle a 
Fenn 1979 

Sea salt 

1.500- 2 x10'® 

1.500- lx 10‘® 

1.490 - 2x 10*^ 

1.470-2x10*® 

1.456- 7x 10*® 

1.470 - 0.0025 

1.500 - 0.014 

Shettle £ 
Fenn 1979 

Crustal 

- 7x lo'^ 

- 6 X 10‘^ 

- 2.5x 10*^ 

- 2x 10*^ 

- 2.5 X 10*’ 

- 2 X lO"^ 

- 0.2 

Patterson 

1981 

75% 

1.432 - 2 X 10*® 

1.431 - 2 X lo'® 

1.429 - 2 X 10*® 

1.420- 1.5x 10*® 

1.398 - 2.72 X 10*® 

1.359 - 0.123 

1.737 - 0.273 

Palmer £ 
Williams 197! 


1.77- 2x lo"^ 

1.77 - 2 X lo"’ 

1.77 - 2 X 10*’ 

1.76- 6x 10*® 

1.74 - 5.5 X 10*® 

1.62 - 3.1 x 10“^ 

0.55 - 0.061 

Toon & 
Pollack 1976 

NaCl 

1.55 - 10*^ 

1.55 - lO"’ 

1.55 - 10*^ 

1.53 - 10*’ 

1.53 - 10*’ 

1.515 - lO"’ 

1.49 - 10“^ 

Toon £ 
Pollack 1976 

(NH^)2S0^ 

1.53 - lO"^ 

1.52 - lO"’ 

1.52 - 10*’ 

1.51-2.4x10*® 

1.49- 1.8x 10*® 

1.46 - 0.006 

1.98 - 0.06 

Toon £ 

Pollack 1976 

Sahara 

dust 






1.56 - 0.015 

1.74 - 0.40 

Volz 1973 

Volcanic 

dust 






1.56 - 0.009 

1.95 - 0.40 

Volz 1973 




Figure 3. 1 Schematic diagram showing the variation of aerosol properties with altitude. 
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Figure 3.2 Latitude distribution of the variation of stratospheric aerosols {from Rosen et 
ai. 1975). Numbers on the diagram give the number of aerosol particles per mg of air. 
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Figure 3.3 Typical samples of aerosol size distribution in different regions of the at- 
mosphere: (a) boundary layer; (b) free troposphere: (c) stratosphere. 
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Figure 3.4 Location of free tropospheric and stratospheric aerosol measurements used in 
this report. 
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PARTICLE RADIUS. 


Figure 3.5 Aerosol volume distributions measured over the USA by Blifford and Ringer 
(1969). 



Figure 3.6 Aerosol volume distributions measured over Europe by Cress (1980). 
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Figure 3. 7 Aerosol volume distribution measurements made over the USA or the Pacific 
Ocean by Patterson et al. (1980). 



Figure 3. 8 Volume size distributions measured by different research workers over land and 
ocean (Blifford, 1970; Savoie, 1978; Patterson et al., 1980). 
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Figure 3.9 Altitude variation of large particle concentrations (r > 0.2 ijun) over Europe 
(Cress, 1980). 
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Figure 3.10 Altitude variation of large particle concentrations (r > 0.5 \jon) over the USA 
(Patterson et al., 1980). 
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Figure 3.11 Altitude variation of large particle concentration (r > 0.5 pm) over the Pacific 
Ocean (Patterson et al., 1980). 



LATITUDE, DEGREES 


Figure 3. 12 Latitude variation of aerosol mass loading in the free troposphere over the 
Pacific Ocean (Patterson et al., 1980) and the equivalent /Jcoj. 
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Figure 3.13 Cumulative probability distribution for the aerosol particle concentrations 
over Europe (Cress, 1980). N/N„ represents the ratio of the aerosol number concentration 
to the mean concentration. 
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Figure 3. 14 Cumulative probability distribution for the aerosol particle concentration over 
the USA and the Pacific Ocean (Patterson et al., 1980). N/N„ represents the ratio of the 
aerosol number concentration to the mean concentration. 
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Figure 3. 15 Cumulative aerosol size distributions obtained by Hofmann and Rosen (1981). 





4. OPTICAL MEASUREMENTS OF AEROSOLS 


4.1 Lidar Measurements 

Lidar measurements have now been made for approximately twenty years and have 
developed from very simple experiments, designed to record total molecular and aerosol 
backscatter, to a very wide range of techniques for studying different atmospheric species 
and parameters. Of particular interest to this report are measurements of aerosol backscat- 
ter. These measurements, in one form or another, go back to the beginning of lidar as a tech- 
nique and a small number of data sets exist that extend over many years. The most consistent 
and reliable of these are those for stratospheric aerosols and the records that exist clearly 
show the aerosol increases due to volcanic eruptions and the subsequent slow decreases. 
Tropospheric data, although extending back as far as stratospheric, presents a much less 
orderly picture. There are two reasons for this. In the lower troposphere, large variations 
occur in aerosol concentration. Relatively few systematic comparisons have been made of 
lidar returns with aerosol source and history and there are no long-term comparisons. In the 
upper troposphere the aerosol concentration is low. Unfortunately, lidar systems are dif- 
ficult to calibrate directly and the returned signals are most commonly normalized by assum- 
ing that the scattering is entirely molecular at a selected altitude in the atmospehre. For 
tropospheric measurements, the altitude most often chosen is that just below the 
tropopause, rendering it impossible to assess the aerosol contribution at that height. 

Early lidar measurements were made using ruby lidar systems operating at a wavelength 
694 nm. Ruby systems are still much used although they are considerably supplemented by 
Neodymium Yag at 1.06 ^m and, much less commonly, by harmonics of those wavelengths. 
A few recent observations exist at 10.6 ijm but virtually no observations have been made for 
wavelengths between 1 .06 and 10.6 /im. The geographic spread of lidar observations is poor, 
being heavily concentrated in North America, Europe, and Japan. This is shown in Fig. 4.1, 
which indicates those stations making tropospheric and lower stratospheric aerosol 
measurements using visible or near-infrared wavelengths. 

In order to assess the situations as regards available data sets on aerosol scattering using 
lidar, letters were written to research workers believed to have made such measurements, 
particularly those outside the USA. Table 4. 1 summarizes the replies and other information 
gained, either by direct contact or from published reports (Names and addresses of those 
persons supplying information are included in the appendix to this report.) Section (a) of the 
table is a very brief summary while Section (b) gives more detailed information. For the sake 
of completeness the table also contains information about CO 2 lidar measurements on 
aerosols and lists stations making measurements on fog, cloud, etc. It may be noted that 
only one research group in the USA, namely, that of Marshall Space Flight Center, is at pre- 
sent making airborne measurements of CO 2 wavelengths. Stations making lidar 
measurements on gaseous atmospheric species only, are not listed. 
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Earlier in this Section it was noted that a few good long-term data sets existed of lidar 
measurements of the stratospheric aerosol. One such, probably the best, obtained at Hamp- 
ton, Virginia, is shown in Fig. 4.2 (Swissler et al., 1982). The figure shows the peak backscat- 
ter cross section at a wavelength of 0.6943 fjm between 1974 and 1979. (The scale on the 
right-hand side of the figure which shows the equivalent backscatter at a wavelength of 10.6 
/um will be discussed in Section 7.) The effects of the eruptions of Fuego and Mount St. 
Helens can be seen very clearly. Table 4.2 shows a list of those volcanoes in the past three 
years known to have ejected material into the stratosphere. This has been a period of very 
considerable volcanic activity, culminating in April 1982 in the eruption of El Chinchon, 
probably the most significant eruption (as far as the stratosphere is concerned) since that of 
Agung in 1963. A short table listing the most important conclusions concerning lidar 
measurements is included in Section 8. 

4.2 Extinction Measurements 

Extraterrestrial radiation passing through the atmosphere will be absorbed and scattered 
by air molecules (including ozone, water vapor, CO 2 and other absorbing species) and at- 
mospheric aerosols, resulting in a reduction in intensity according to the Lambert-Beers law 
(at a single wavelength). Since the extinction in radiation is directly related to the optical 
properties of the aerosol, the aerosol size distribution and/or its concentration can be ob- 
tained by inverting the measured extinction, provided we can assume some other properties 
of the aerosol, such as its index of refraction. 

Atmospheric aerosol extinction measurements can be divided into two general classes: in 
the first the total extinction due to all aerosols along the path between the sources and the 
detector device is measured; in the second the extinction from aerosols along a path of unit 
length is measured or deduced. Atmospheric turbidity, visibility and total optical depth 
measurements belong to the first class. These measurements are relatively easy to make and 
so provide good continuity. However, results of these measurements are usually related to 
the mean aerosol properties, and the distribution of aerosol properties along the path of 
radiation is almost impossible to obtain from the experimental data. On the other hand, 
those measurements belonging to the second class, such as extinction or attenuation coeffi- 
cient measurements, conducted on the ground or by detection devices mounted on space 
craft can provide information on aerosol properties with, in the latter case, a resolution of 
about 1 km. For the purpose of this study we are more interested in the extinction data ob- 
tained by measurements of this second class mentioned above. 

Although aerosol extinction or attenuation coefficient profiles obtained by numerical 
modeling have been reported in the litrature (Elterman, 1970; Kneizys et al., 1980), we are 
aware of only a very few measured extinction profiles (e.g., Elterman, 1966). So far the 
largest data base on aerosol extinction coefficients is undoubtedly that obtained by the SAM 
II and SAGE satellite experiments. The SAM II instrument consists of a single channel sun 
photometer, centered at 1.0 ^im wavelength, mounted on the Nimbus-7 satellite which was 
launched on October 23, 1978. Aerosol extinction coefficient profiles, with a 1 km vertical 
resolution, are being obtained by this method. The SAGE instrument consists of four 
radiometric channels centered at 1.0 mu, 0.6 /urn, 0.45 >im, and 0.385 ^m, mounted on the 
Application Explorer Mission-B (AEM-B) satellite, which was launched on February 18, 
1979. The four-channel extinction measurement obtained by SAGE can be inverted to ob- 
tain aerosol extinction coefficients at 1 .0 /.<m and 0.45 /urn. Due to the orbital characteristics 
of the Nimbus-7 satellite, SAM II measurements are limited to two latitude bands ranging 
from 64°S to 81°S, and 65°N to 83°N. In contrast, the highly processing orbid of the AEM-B 
satellite allows the SAGE measurements to cover global areas in the latitudes between about 
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79°S and 79°N. Though the majority of data obtained by both SAGE and SAM II is con- 
fined to stratospheric regions below 30 km, under favorable conditions without cloud 
coverage, extinction data to the surface of the earth can sometimes be inverted. If the con- 
version of SAGE and SAM II extinction data to Pco 2 is proven to be feasible, the spatial and 
temporal variation of Pco^ over the globe can be accurately modeled. 
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TABLE 4.1: LIDAR MEASUREMENTS ON AEROSOLS, FOG AND CLOUD IN THE LOWER ATMOSPHERE 


(a) Summary 



(i) 

MADE USING VISIBLE OR NEAR INFRARED WAVELENGTHS 
Aerosol Measurements 

(ii) Fog, Cloud or 
Precipitation Only 

MADE USING 
CO^ WAVELENGTHS 



Available Data Sets 







No. of 

No. of 

Extensive 

Limited 

Extensive 

Limited 

No. 

of stations 

No. of No. of 

Continent 

stations 

airborne 

tropospheric 

tropo- 

stratospheric 

stratos- 



stations airborne 





spheric 


spheric 



systems 

North America 

13 

2 

4 

5 

1 

3 


4 

3 1 

and Pacific 










South America 

2 





2 




and Caribbean 










Europe 

8 

2 

1 

7 

1 

2 


3 

4 1 

Asia 

6 


1 

4 

1 



2 


Australia 

2 


1 

1 

1 






Notes: 1) Stations and airborne systems no longer operating are not listed unless there exists a 

considerable data base. 

Extensive stratospheric observations imply a data base of several years; extensive 
tropospheric observations imply several sets of experiments. 


2 ) 








(b) Detailed Information 


Station Principal Wavelengths Period of Source of Information Notes 

Research Used Operation 

Worker (s) 


North America 
and Pacific 



U. of Wisconsin 

E. 

Eloranta 

694 nm 


Conference Proceedings , 

Extensive studies of convec- 


Meteorol. Dept. 

J. 

Weinman 



Publications 

tive boundary layer using 


Madison, WI 





Phone conversation 

scanning lidar system. High 







Nov. 30, 1981 

spectral resolution airborne 
lidar for separation of 
aerosol and molecular 
scattering. 


Air Force Geophysics 

D. 

E. Bedo 

1,064,355 nm 


Conference Proceedings 

Proposed balloonborne lidar 


Laboratory, Hanscom 

K. S.W. Champion 



system for high altitude 


AFB, MA 01731 

U.S.A. 






studies. 


Photometries Inc. 

G. 

Davidson 

0.53 ym 

1981 

Conference Proceedings 

Backscatter measurements in 


Woburn, MA 01801 






falling snow. 


Atmos. Sci. Center 

E. 

Uthe 

694 nm 

1960s 

Publications, reports 

1) Numerous lidar systems 


SRI International 

P. 

Russell 

(also CO^ 

present 


developed including mobile 


Menlo Park, CA 



laser) 



and airborne 


94025, USA 2) Main field- of study is 

tropospheric, particularly 
boundary layer 

3) Shipborne measurements 
made in 1977 in Atlantic 
and Mediterranean 

4) Considerable work done on 
modelling of aerosols. 



Station 

Principal 
Research 
Worker (s) 

Wavelengths 

Used 

Period of 
Operation 

Source of Information 

Notes 

NASA- I^ang ley 
Research Center 
Hampton , VA 
23665, USA 

E. V. Browell 

600 nm 

1980 

Conference reports, 
personal contact 

Airborne system has been used 
to obtain downward looking 
profiles from altitudes of 
about 3 km. Aerosols, ozone, 
and atmospheric water vapor 
studied. 

NASA-MSFC 

Huntsville, AL 35802 

J. Bilbro 

10.6 pm 

1981-82 

Conference proceedings, 
personal contact 

Airborne aerosol backscatter 
and doppler shift measurements. 

NASA/GoddarcJ 
Space Flight Center, 
Code 913, 

Greenbelt, MD 20771 

J. D. Spin- 
hirne 

532 nm 

1979-82 

Conference proceedings 

Investigation of cloud 
structure . 

U. of Arizona 
^ Dept, of Elec. Eng.* 

or 

Institute of Atmos. + 
Physics 
U. of Arizona 
Tucson, AZ 85721 
U.S.A. 

*J. Reagan 
+R. Schotland 

694 nm 
8-12 ym 

1970 -*■ 
present 

Publications, 
personal contact 

Monstatic and bistatic lidar, 
several experiments conducted 
in and closely above the 
boundary layer combining 
several observational 
techniques. Recent CO 2 lidar 
measurements . 

Atmospheric 
Sciences Lab 
White Sands Missile 
Range, NM 88002 

J.S. Randhawa 
E. Measure 

694 nm 
(also CO 
lidar) 

1979 -»■ 
present 

Conference reports, 
personal visit 

Application mostly to 
smoke and dust, ongoing 
data collection comparing 
visible and CO^ backscattering 
and extinction. Measurements 
in U.S.A. and Europe. 

U.S. Army Atmospheric W.J. Lentz 

Sciences Laboratory 
Atmospheric Sensing 
Division 

White Sands Missile Range, 

NM 88002 

1.06 ym 

1981-82 

Conference proceedings 

Main interest is in inversion 
of data to obtain extinction 
profiles. 


I 



station Principal Wavelengths Period of 

Research Used Operation 

Worker (s) 


Jet Propulsion R. T. Menzies 10.6 ym 

Laboratory 

California Institute of 
Technology 
Pasadena, CA 91109 

Environmental Monitor- J.L. McElroy 1.06, 0.53 Urn 
ing Systems Laboratory 
U.S. Environmental 
Protection Agency 
P. 0. Box 15027 
Las Vegas, NV 89114 

NASA-Langley M,P. McCormick 694 nm 1974 - 

^ Research Center T.J. Swissler (ground-based) present 

-ii Hampton, VA 23665 J. Goad 694 & 1060 nm 

(airborne) 


Environmental V. Derr 694 & 347 nm 1973 - 

Research Labs. J. deLuisi (also CO present 

NOAA F. Hall lidar) 

Boulder, CO 80303 W. Eberhard 

R. Schweisow 


Aeronomy Laboratory C.S. Gardner 596 nm 1980 

Dept, of Elec. Eng. C.F. Sechrist 

U. of Illinois 
Urbana, IL 61801 

Barrow, Alaska J. deLuisi 

c/o NOAA 

Env. Res. Labs. 

Boulder, CO 80303 


Source of Information 


Notes 


Conference proceedings, Lidar returns from aerosols 

personal contact and topographic targets. 


Conference proceedings Downward looking airborne 

lidar to study air pollution 
transport . 


Numerous publications 1) Good back- log of data on 

and reports, personal stratospheric aerosol 
contact 2) Detailed study of 

stratospheric volcanic effects 

3) Operational airborne 
system, upward & downward 
looking from aircraft 

4) Polarization studies. 

Stratospheric data following 
Mt. St. Helens eruption. May 
1980 (CO^ and visible lidar 
comparison) . Studies of plume 
dispersion CO aerosol measure- 
ments from 1980, including 
backscattering and doppler shif- 

Studies of mesopheric 
sodium, stratospheric aerosol 
observations following Mt. St. 
Helens eruption. May 1980. 

Recent correspondence System set up but not yet 
(November 1981) operational, long-term 

future uncertain. 


Conference presen- 
tation; recent corres- 
pondence (Dec. 1981) 


Publications and 
conference presen- 
tation, personal 
contact 



station 


Principal 
Research 
Worker (s) 


Wavelengths 

Used 


Period of Source of Information 
Operation 


Notes 


NOAA Air Res. 

Laboratory 
Hilo, HI 95616 

J. deLuisi 
R.W. Pegley 
F, Fernald 

694 nm 

1973 - 
present 

Recent correspondence 
(November 1981) 
and conference 
presentation 

Several years data, processing 
incomplete; continuing to 
monitor stratospheric aerosols. 

University of Utah K. Sassen 

Dept, of Meteorology 
Salt Lake City, UT 94112 

694 nm 

1982 

Conference proceedings 

Studies of scattering 
characteristics of clouds 
(including polarization) . 

Defense Research R.E. Kluchert 

Establishment, 

Valcentier 
P. 0. Box 8800 
Courcelettev Quebec , Canada 

1.06 lim 


Conference proceedings 

System used to map 
pollution plumes 

York University 
4700 Keele St. 
Pownsview, Ontario, 
Canada M3J 1P3 

A. Carswell 


1977-82 

Publications and 
conference proceedings 

Studies of scattering in 
clouds and fogs. 

Defense Research 
Establishment, 
Valcentier 
P. 0. Box 8800 
Courcelette, Quebec, 

E . T . Evans 
Canada 

1.06 pm 


Conference proceedings 

Lidar mapping a 
of smoke aerosols. 

South Rraeilca 
and Caribbean 






INPE-C.P. 

Sao Jose dos Compos 

Sao Paulo 

Brazil 

B.R. Clemesha 

694 and 
589 nm 

1970 - 
present 

Recent correspondence 
(November 1981) 

Continuous data base to 1977 on 
stratospheric aerosols, 
B-profiles for 1970-77 
between 10 & 50 km published. 

Dept, of Physics 
U. of the West Indies 
Kingston 7 
Jamaica 

M. Shams 

694 nm 

1965-1980 

Personal operation 

Intermittent data on strato- 
pheric aerosol for 1965- 
1980. Total data base is 
extensive. Present status 
uncertain. 



station 


Principal 
Research 
Worker ( s ) 


Wavelengths 

Used 


Period of 
Operation 


Europe 

Instituto di Ricera L. Stefanutti 1.06 um 
sulle Onde Electro- P. Bruscaglioni 0.53 pm 
magnetische of CNR 
Via Paniatichi 64 
Firenze 50127 
Italy 

National Swedish K. Fredriksson 1.06 pm and 

Environment dye laser 

Protection Board 
S-611 82 Nykdping 
Sweden 

Department of Elect- S. Lundquist 10.6 pm 

rical Measurements 
Chalmers University of 
Technology 
S-412, 96 GUteberg 
Sweden 

National Defence A. H3gard 9.2 ->■ 10.8 pm 1982 

Research Institute 
Department 3 
Box 1165, S-581 11, 

Linktthing, Sweden 

University of Bergen I. Singstad 694 nm 1982 

Norway 


1975 - 
present 


1978 - 
present 


Physics Lab TNO 

The Hague 

The Netherlands 


C.W. Lamberts 1.06 pm 


1978 - 
present 


Source of Information 


Notes 


Conference reports; Main research has been on 
recent correspondence properties of fogs. Data 
(Dec, 1981) base on aerosols limited. 


Conference reports; Main interest is in urban 

recent correspondence pollution (so , NO^, 0 , HO). 
(Dec. 1981) No absolute measurements of 

aerosol cross-section. 


Recent correspondence Measuring diffuse 

(Feb. , 1982) reflectivity coefficients. 


Conference proceedings. Aerosol scattering and 

personal contact extinction from fog, snow, 

rain, etc. 


Conference proceedings Measurement of Stokes Parameters 

in scattering from clouds. 

Conference reports; Program of measurement of 

recent correspondence extinction and backscatter 

(Dec. 1981) caused by tropospheric 

aerosols. Future work comparinc 
1.06 pm and 10.6 pm scattering 
and extinction planned. 
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Station 

Principal 
Research 
Worker (s) 

Wavelengths 

Used 

Period of 
Operation 

Source of Information 

Notes 

Service d'Aeronomie 
BP3 91370 
Verrieres-Buisson 
France 

A. Hauchecorne 
L. Chanin 
G. Megie 

1.06 pm, 
530, 353 nm 
and various 
dye lasers 

1970 - 
present 

Publications and 
conference reports; 
recent correspondence 
(Feb. 1982) 

An extensive range of 
measurements is being 
conducted on a wide range 
of atmospheric properties. 
Information on scattering 
from volcanically enhanced 
stratospheric layers has 
been published but no other 
aerosol data base. 

Direction de la 
Meteorologie 
-EERM 

78470 MAGNV Les Hameaux 
France 

J.L. Gausset 

t 

694 nm 

1982 

Conference proceedings 

Measurements of visibility 
in fogs. 

Meteorologie D. Renaut 

Nationale, EERM, 

78470 St-Remy-les-Chevreuse 
France 

347,266 nm 

1981 - 
present 

Conference proceedings 

Raman measurements of 
aerosols, water vapor and 
gaseous pollutants in the 
boundary layer. 

Deutscher Wetterdienst 
Meteorologische 
Observatorium Hamburg 
Frahmredder 95 
2000 Hamburg 65 
FRG 

F. Kasten 
E . Klaphack 


1979 - 
present 

Conference report, 
letter written; 
recent correspondence 
(Nov. 1981) 

Studies of tropospheric 
and anthropogenic aerosols 
(0-2 km altitude) using 
steerable lidar system. Data 
not yet available. 

DFVLR - Inst, for 
Optoelectronics 
D-8031 Oberpfaffenhofen 
FRG 

Ch. Werner 
F. KfJpp 

1.06 ym 
(airborne 
system) 

694 nm 

ground based, 

10.6 pm 

1970 

Conference reports, 
publications, personal 
experience; recent 
discussions (June, 1982) 

Extensive observations and 
system development over a 
period of many years. Interest 
mainly on tropospheric 
aerosols. Operational airborne 
system. CO system 
in use for doppler shift 
measurements. Ship-borne 
studies of plume dispersion. 
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Station 


Principal 
Research 
Worker (s) 


Wavelengths 

Used 


Period of 
Operation 


DFVLR - Inst, for W. Renger 

Atmospheric Optics 
D-8031 Oberpfaf fenhofen 
FRG 


1.06 pm, 
530 nm in 
airborne 
system 


Fraunhofer Institute R. Reiter 694, 347, 530 1976 - 

for Atmos. Env. Res. nm present 

D-8100 Garmisch- 
Partenkirchen, FRG 


SRC Rutherford & L. Thomas* 694 nm 1966 - 

Appleton Laboratories and dye present 

Ditton Park lasers 

Slough, Berks SL3 9JX 
England 

Royal Signals and J. Vaughan 10.6 ym 1981-82 

Radar Establishment, 

St. Andrews Road 
Great Halverson, 

Wores WR13 3PS 
England 

Department of E.L. Thomas 10.6 ym 

Applied Physics, 

University of Hull 
Hull, 

United Kingdom 


* Now at The University College of Wales, Aberystwyth, Wales. 



Conference reports. Operational downward looking 

personal discussion; airborne system for the study 

recent correspondence of aerosols and pollution. 

(March, 1982) Measurements over Alpine 

regions . 

Published data on back- 
scattering profiles for both 
tropospheric (1-3 km) and 
stratospheric (10-30 Ion) 
aerosols. Study of volcanic 
effects. Making extensive, 
multi-wavelength transmission 
measurements in the boundary 
layer; data not yet published. 

Series of experiments over 
many years. Main interest 
stratosphere and mesosphere . 
Some data on aerosol 
scattering published. 

Recent correspondence Considerable body of back- 

(Feb., 1982) scatter data from airborne 

system. Not yet published. 


No absolute measurements of 
volume backscattering 
coefficient yet, but future 
measurements planned. 


Conference proceedings ; 
publications; recent 
correspondence (Nov. 
1981) 


Publications; 
recent correspondence 
(March, 1982) 


Recent correspondence 
(Nov., 1981) 
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station 


Principal 
Research 
Worker (s) 


Wavelengths 

Used 


Asia and Australia 

The Institute of V.E. Zuev 1.06 ym 

Atmospheric Optics 530, 490 nm 

U.S.S.R. Academy of 
Sciences, Siberian 
Branch, Tomsk 
U.S.S.R. 

Institute of X.J. Zhou 694 nm 

Atmospheric Physics Y.Z. Zhao 

Academic Sciences 
Beijing, 

China 


Water Research Inst. 

K. Nagaya 

694 nm and 

Nagoya University 

Y. Iwasaka 

other 

Chikusa-hu 


wavelength 

Nagoya 



Japan 




Dept, of Physics M. Hirono 1.06 pm 

Kyushu University M. Fujiwara 530 nm 

Fukuoka 812 
Japan 


Research Institute H. Inaba 

of Elec. Communications 
Tohuku University 
Katshiva 2-1-1 
Sendai 980 
Japan 


Period of Source of Information 
(iteration 


tfotes 


1972 - Conference reports; 

present recent correspondence 

{ April, 1982 ) 


1972 - Conference reports; 

present recent discussion 

(June, 1982) 


1979 - Conference reports; 

present recent correspondence 

(March 1982) 


1972 - Publications and 

present conference presen- 

tations; recent 
correspondence (March, 
1982) 

Conference proceedings 


Extensive work carried out on 
scattering and extinction by 
tropospheric aerosols. 


Interest is in tropospheric 
aerosols and their extinction. 
Extent of data base unknown. 
DIAL measurements of atmospheric 
water vapor. Recent visibility 
measurements . 

Interest in tropospheric 
aerosols 0-2 km. Published 
data on aerosols from dust- 
storms over China. Strato- 
spheric data available also. 

Has good continuous data 
base on stratospheric aerosols. 
Planning measurements of 
aerosol profiles from 6 km 
upwards • 

Main interest is in water 
droplets in fogs and clouds. 
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Station 


Principal 
Research 
Worker (s) 


Wavelengths 

Used 


Period of Source of Information 
Operation 


Notes 


National Inst, for 
Environmental Studies 
Yatabe, Tsukuba 
Ibaraki, 305 
Japan 

N. Takeuchi 

1.06 m 
532 nm 

1979? - 
present 

Conference presen- 
tations; recent 
correspondence (March, 
1982) 

Main interest is in spatial 
and dynamic structure of air 
pollution. Using large 
atmospheric multi-purpose lidar 

Department of Atmos- 
pheric Sciences 
The Hebrew University 
of Jerusalem, 
Jerusalem, Israel 

A. Cohen 



Publications, 
conference proceedings 

Studies of scattering in the 
lower atmosphere and in 
clouds . 

Space Physics Division 
Vikram Sarabhlai space 
Centre, 

Trwinitrum - 695022, 
India 

B.V. Krishna 
Mur thy 

514. 5 nm 

1982 

Conference proceedings 

Bistatic aerosol studies. 

Dept, of Physics 
U. of Adelaide 
Adelaide 
Australia 

S.A. Young 
W.G. Elford 

694 nm 

1969-73 
(at least) 

Publications , 
personal contact; 
recent correspondence 
(Nov. , 1981) 

Good data base on scattering 
from tropospheric and 
stratospheric aerosol. No 
longer operational. 

CSIRO Division 
of Atmos. Physics 
Aspendale 
Victoria 
Australia 

C.M.P. Platt 

694 nm 

1971 - 
present 

Publications , 
conference and 
workshop reports , 
personal contact; 
recent correspondence 
(Feb., 1982) 

Main interest has been on 
scattering properties of 
clouds. Some boxindary layer 
backscatter cross-sections. 
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TABLE 4.2: RECENT VOLCANIC 


DATE 

VOLCANO 

April 17. 1979 

Soufriere 

November 13. 1979 

Sierra Negra 

May 13. 1980 

St. Helens 

Oct 7. 1980 

Ulawun 

April 27. 1981 

Alaid 

May 15. 1981 

Pagan 

Dec-Jan. 1982 

Possibly Pagan 

March 19. 1982 

St. Helens 

March 28 - 
April 4. 1982 

El Chinchon 


AFFECTING THE STRATOSPHERE 


LOCATION 


COLUMN HEIGHT 

13.3 

N. 

61.2 W 

18 - 

20 

km 

CO 

o 

S. 

91.2 W 

> 

14 

km 

46.2 

N. 

122.2 

W 

> 

23 

km 

5.0 

S. 

151.3 

E 


10 

km 

50.8 

N. 

155.5 

E 


15 

km 

18.1 

N. 

145.8 

E 

13 - 

20 

km 


Lidar at 10 N 
shows peak at 17 km 

122. 2W 14 km 


17.3 N. 93.2 V! 


26 km 




Figure 4.1 Map showing lidar stations making tropospheric and lower stratospheric 
aerosol measurements using visible or near infrared wavelengths. Lidar stations are marked 
bya -tt. 



TIME 


Figure 4.2 Stratospheric lidar record at Hampton, Virginia (Swissler et al., 1982) and the 
equivalent 10.6 yon backscatter. 
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5. TECHNIQUES FOR THE CALCULATION OF Pco, 

5.1 Introduction 

It is possible to model Pco^ by the following three techniques: 

(a) Direct calculation of the /Jcoj values from the reported aerosol properties including 
size distribution, concentration, and composition or complex refractive index. 

(b) Extrapolating the existing backscattering data at visible or near infrared wavelengths 
to CO 2 laser wavelengths (e.g., A = 10.6 pm). 

(c) Converting the existing aerosol extinction data to backscattering values at the same 
wavelength or, more directly, to the CO 2 wavelength that we are most interested in. A more 
detailed discussion of these three methods will be given in subsections 5.2, 5.3, and 5.4, 
respectively. 

In our calculations, we assume that the aerosol particles are spherical so that a Mie code 
can be used. The assumnption is justified because most of the reported aerosol particles in 
the atmosphere have size ranges from about 0.01 pm to a few microns. Particles of that size 
range are usually spherical with the possible exception of those larger particles close to a few 
micron radius. The errors involved by assuming the sphericity of these larger particles will be 
discussed in Section 6.2.4. 

The majority of our calculations have been carried out by the scattering code AGAUS 
developed by the Atmospheric Sciences Laboratory of US Army at White Sands Missile 
Range, New Mexico. As a check, the results of this code have been compared with that of an 
older scattering code developed by N. Logan. We are completely satisfied with the AGAUS 
code because results from their code are different from Logan’s results by less than about 
1%, even when integrations over complete size distributions are carried out. The original 
AGAUS code can only accept commonly used analytical expressions for aerosol size 
distributions or else numerical concentrations specified at equal size ranges. It outputs 
several optical properties, including extinction coefficients and backscattering functions. In 
order to suit our purpose we have slightly modified the original program so that it can accept 
measured aerosol size distributions in formats reported in open literature as well as use dif- 
ferent methods to extrapolate the input data to larger sizes. 

Since large variations occur in aerosol size distributions, most of our analyses have been 
conducted with three approaches. The first is to conduct calculations using a log-normal size 
distribution and different mode radii. A semi-quantitative conclusion can be drawn from 
this aproach. Our second approach has been to use some typical aerosol size distributions 
measured in different regions of the atmosphere. The third is to calculate the backscatter 
functions for a large number of possible size distributions and interpret the values obtained. 

5.2 From Measured Aerosol Size Distributions and Models 

The two input parametrs of ftco 2 calculation are the size distribution (including total con- 
centration) and the refractive index of the aerosol material at the CO 2 wavelength. 
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As mentioned earlier, our first approach is to assume aerosol particles have a log-normal 
size distirbution. The width of the mode is fixed at 2.0 /urn and the mode radii range from 
0.05 ^m to 5.0 /urn. Our second approach is to use the aerosol size distributions obtained by 
in situ measurements. A list of some of the reported aerosol size distributions used in this 
study is given in Table 5.1. From the size range listed in this table, it can be seen that, with 
the exception of Cress’ data, most measurements are cut off at a few microns particle radius. 
In order to study the effect of large particles on the value of ftcoj, we have also calculated 
Pco 2 with aerosol size distributions extrapolated to 10 /im, or even larger sizes, using both 
linear and log-normal extrapolation. A comparison of the resultant /Jcoj values will be 
presented in the next section. 

Some of the reported aerosol size distributions are not smooth but show fluctuations in 
the larger size range. In some of these cases we have drawn envelopes to the reported aerosol 
size distribution and calculated the corresponding lower and upper limits of /Icoj- 

The aerosol materials under investigtion are sulfuric acid solution with 75Wo H2SO4 by 
weight, water soluble, dust, soot, clay, and (NH4)2S04. Their corresponding indices of 
refraction at 10.6 j^m are listed previously in Table 3.5. 

5.3 From Lidar Measurements 

Lidar measurements made at visible or near-infrared wavelengths may be used to calculate 
the equivalent value for ftco 2 provided a suitable conversion factor for the ratio or 

n.or in/rar.d) Hiay bc fouttd. This factor will depend upon the wavelengths used but, more 
significantly, it will depend upon the refractive indices and size distributions of the aerosols 
under study. These will normally not be known well, and the conversion factor will not be 
exact. The question of whether visible lidar measurements form a useful data base for con- 
version to CO2 wavelengths depends therefore on whether useful limits may be placed on 
this factor whose values vary from 10"* to 10 ', depending upon particle size, wavelength and 
refractive index. In particular, it is relevant to ask whether more accurate results may be ob- 
tained by this technique than by carrying out calculations based directly upon the estimated 
aerosol parameters themselves. This question will be directly addressed in Section 6.2.2. 

5.4 From Extinction Measurements 

The conversion of extinction measurements, made at GO2 or other wavelengths, to 
equivalent backscattering function at CO2 wavelengths is essentially similar to the conver- 
sion of lidar measurements made at other wavelengths. The caveats with regard to the 
establishment of a suitable conversion factor, discussed in relation to lidar in the previous 
section, apply equally well to extinction. Our conclusions with regard to this problem will be 
discussed in Section 6.2.3. 

The problem with regard to the use of optical depth measurements is more complex. Op- 
tical depth is equivalent to the total extinction integrated between the two points of 
reference. For many measurements this is the whole atmosphere from space to the point of 
measurement. The extinction coefficient will vary along the path and the characteristics of 
aerosol, as shown by the optical depth measurement, will be heavily biased toward the 
region of maximum extinction. For example, multi-wavelength radiometer measurements 
made at the earth’s surface are used to deduce the mean aerosol particle size distribution 
along the path (Reagan et al., 1977). Since, in many cases, the maximum aerosol concentra- 
tion will be near the point of measurement, the size distribution found will correspond most 
closely to this region, and give relatively little information about the aerosol at other points 
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along the path. This is not to say that optical depth measurements may not give useful infor- 
mation about the aerosol in the free troposphere when measurements are made from a point 
about the boundary layer (e.g., Shaw, 1980). Once the region in which the extinction is oc- 
curring has been established, conversion to Pcoi be attempted, with additional inac- 
curacies arising because it is an optical depth rather than an extinction measurement. This 
necessitates assumptions about the size of the region and the variations of relative extinction 
within the region. 
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TABLE 5 . 1 ; Examples of reported 
used in this study 


Date 


Region 
and Height 


Type 



Laramie , 
Wyoming 
16 km 


Mildura, 
Australia 
10, 16, 23 km 

Sherborne , 
England 
1-6 km 

U.S.A & 
Pacific 
0-6 km 


stratospheric 


stratospheric 


European 

continental 


Size Range 


Author 


0.01 ym to 1 ym, Hofmann and Rosen (1981) 

but larger size 
aerosols exist 

0.05 ym to about Gras and Ayers (1979) 

0.9 ym 

~0.2 ym to ~8 ym Cress (1976) 


0.1 to 3 ym 


Patterson et al (1980) 













6. ERRORS AND UNCERTAINTIES IN THE CALCULATION OF Pco, 


6. 1 Sources of Error 

Since it is not feasible to obtain complete and accurate information about all the relevant 
properties of aerosols for /Jcoj modeling, errors are expected to occur. It is therefore essential 
to identify the sources of error and estimate the magnitude of the uncertainties incurred 
from such errors. 

We can identify three types of errors; those associated with the measuring instruments, 
those associated with the analysis technique, and those associated with the properties of the 
aerosol particles themselves. The possible errors associated with the measuring instruments 
have been discussed in Section 3.2. Here we just wish to point out that, due to the instability 
of the measuring systems, noise may occur in the measured aerosol property, especially in 
the size distribution. As also pointed out in that section, some in situ measurements are ef- 
fectively not in situ; the properties of aerosol particles that we are interested in have to be 
deduced in the laboratory from the samples collected in the field. As will be discussed later in 
this section, different aerosol properties can be deduced from the same measured results 
using different analyzing techniques. The incompleteness and inaccuracy of the reported 
aerosol properties are another source of error. Due to the limitation of conventional measur- 
ing devices, the small number concentration of particles in the coarse-particle mode is hard 
to measure. Consequently, their properties are usually not reported along with those of 
smaller sizes. This is not a serious problem for modeling extinction and backscatter at visible 
wavelengths, due to their expected small contribution to the overall value. However, for 
radiation of longer wavelengths (such as CO 2 laser), the presence of larger particles may 
greatly increase the overall Pcoi values, as will be reported later in this section. The in- 
completeness of aerosol size distributions can be partially solved by extrapolating the 
reported size range to larger sizes if that seems feasible. However, error may be generated by 
using a wrong extrapolation method. The composition (or more directly the complex refrac- 
tive index) and shape are other sources of error that will also be discussed later in this sec- 
tion. 

6.2 Effects of Errors on the Calculated pco: 

As explained in Section 5.1, three main techniques have been used to estimate errors and 
to assess the dependence of ftco 2 upon an exact knowledge of the parameters involved in its 
calculation. 

The first of these involve the use of log-normal models and the principle features of these 
are summarized in Fig. 6. 1 which shows the key parameters and the relationship between the 
number distribution and the volume distribution. In the following subsections, modeling has 
been carried out for values of mode radius r, between 0.03 jum and 5 pim, covering that of the 
commonly observed range of particle size distributions in the lowest 20 km of the at- 
mosphere. A limited amount of calculation has been done varying the geometric mean stan- 
dard deviation a, between 1.5 and 2.5. 
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Three size distributions have been selected as typical for the free troposphere. These are 
shown in Fig. 6.2. They consist of a continental distrribution with a large number of the 
larger particles (Cress, 1980, spring average), a continental distribution with fewer large par- 
ticles (Patterson et al., 1980) and a marine free tropospheric distribution with very few large 
particles (Patterson et al., 1980). Most aerosol size distributions found in the literature lie 
between the limits represented by these three. 

6.2.1 pcoi calculated from measured size distributions and refractive index 

Using the three size distributions shown in Fig. 6.2 as a starting point, histograms have 
been calculated showing the relative contribution to /Jcoj from different particle size ranges. 
These are shown in Figs. 6.3, 6.4, and 6.5. (Note: the size ranges used in the three figures are 
not the same.) Also shown in these figures are similar histograms for the backscattering 
function at 0.69 and 1.06 ^m and the extinction coefficient at 1.0 fim. The main points 
shown by these figures are: 

(a) That, even for the Cress spring average aerosol size distribution which contains a 
relatively large number of the larger particles, most of the scattering at 10.6 iaiti is from par- 
ticles with radii less than 2.5 >im. In this example, very little scattering comes from particles 
with radii less than 1.0 ^m. 

(b) That, for the marine free troposphere, an appreciable fraction of the scattering comes 
from particles with radii less than 0.5 fim. 

(c) In all cases, the size ranges corresponding to peak scattering at the shorter wavelengths 
differ from those responsible for the peak scattering at 10.6 ixm. 

These histograms and the points noted above will be useful in understanding the source of 
the errors and uncertainties in modeling /Jcoj, to be discussed in the following subsections. 

6. 2. 1.1 Errors related to refractive index uncertainties 

Calculations have been made of the /Jcoj at 10.6 ^m for each of the three aerosol size 
distributions described in Section 6.2 and a range of refractive indices. These are shown in 
Figs. 6.6, 6.7, and 6.8 where the real part of the refractive index is allowed to vary between 

1 . 1 and 2.2 and the imaginary part from 0 to 1 .0. (The figures correspond to the Cress spring 
average, the Patterson et al. continental, and the Patterson et al. maritime aerosol size 
distributions, respectively.) The curves (for constant imaginary component of the refractive 
index) are similar in all three figures, although more overlapping occurs in the Cress spring 
average size distribution than with the other two and the mean level is very different between 
the figures. Also shown on these figures are the scattering and refractive indices for five 
common arosol materials and for water. The following are the main characteristics common 
to these three figures. 

(a) For a zero imaginary component of refractive index /Jcoj varies by about two orders of 
magnitude as the real component of refractive index varies from 1.1 to 2.2. 

(b) For an imaginary component of refractive index of the order of unity (this may occur 
near an absorption band), the scattering is independent of the real component of the refrac- 
tive index. What is in fact happening is that the particle has become totally reflecting. 

(c) For imaginary components of refractive index of less than 0. 1 , the scattering from the 
two Patterson et al. distributions is effectively the same as for a totally nonabsorbing parti- 
cle. 
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(d) If we except water and carbon from the set of aerosol materials shown, the remaining 
four materials have fairly close backscattering functions. Carbon is normally a minor consti- 
tuent although it may be important at certain times and locations. Water is of course not a 
minor constituent and its omission implies that we are dealing with low humidity conditions. 
The importance of the low backscattering functions for water will be taken up again in Sec- 
tion 7.4. The fact that the backscattering functions of the remaining four aerosol materials 
are fairly close implies that errors in refractive index may not have too great an effect on 
/?co 2 . a point that will be discussed again in the next section. 

6. 2. 1.2 Errors related to size distribution uncertainties 

The effect of changing aerosol radius on the values of ftcoi has been examined using log- 
normal models with constant o, = 2.00 and varying r, from 0.03 to 5 ptm. The resultant 
behavior of /Jcoj is shown in Fig. 6.9 for the five different aerosol materials used in Section 
6.2. 1.1. Also shown in Fig. 6.9 are the ranges of log-normal mode maxima found to exist in 
the free troposphere and stratosphere. These have been determined both by actual fitting to 
published size distributions and by noting values quoted by other researech workers. (The 
ranges for the stratospheric and tropospheric accumulation modes are not quite the same 
but have been assumed to be so for the purposes of this figure.) It can be seen that the range 
of Pco 2 is at least six orders of magnitude. Calculation, using actual measured size distribu- 
tions shows a range about as great as this, even though almost all tropospheric aerosol 
distributions consist of a mixture of the two dominant modes. In contrast, as noted in the 
previous section, the change in /Jco^ for different aerosol materials is much less than the 
change produced even by a small shift in mode radius. This point is further emphasized in 
Fig. 6.10 where the values of ftcoi for sulfuric acid, water soluble aerosol (Shettle and Fenn, 
1979) and ammonium sulfate have been plotted as a ratio to ^coj for the dust aerosol (Shettle 
and Fenn, 1979). Since the soluble content of the tropospheric aerosol varies with aerosol 
radius, curves have been plotted for both the four materials and also for a standard mixture 
(Winkler, 1974; Nilsson, 1979). Over the common range of mode radius (.03 1.0 ^m) the 

variation of /Jcoj for the mixed aerosol does not differ by more than a factor of 2'/2 from that 
for a pure dust aerosol, even for the extreme case of ammonium sulphate. In general, 
therefore, it may be concluded that the effects of uncertainties in aerosol composition are 
small compared to those caused by uncertainties in particle size. 

Fig. 6.10 shows two scales, one for the mode radius of the fitted number log-normal 
distribution and one for the equivalent log-normal volume distribution. It may be noted that 
these differ by approximately a factor of four. It is possible to extend this concept and plot 
Pco 2 t not for one particle per unit volume but for unit volume of aerosol per unit volume of 
air. This has been done in Fig. 6.11, and it may be seen that the change in pco 2 is very much 
less than in Fig. 6.9. This is an important result, as it means that an uncertainty in the max- 
imum of a volume distribution will produce less variation in the modelled Pco 2 than a simular 
uncertainty in the maximum of a number distribution. It also means that aerosol volumes or 
mass concentrations may be converted to Pcoi with not more than two orders of magnitude 
possible error, even if no size information is available. There may be a considerable im- 
provement over this figure if a reasonable guess at the type of size distribution is possible. A 
further development of this point is illustrated in Table 6.1. In the calculations described 
above, we have assumed that a, - 2.0. In practice a, for the best fit log-normal distribution 
may vary from 1.5 to 2.5. Table 6.1 shows the effect of this variation in a, on the modelled 
values for /Jcoj- For calculations based on a fixed number distribution mean, the effect is 
severe; for calculations based on a fixed volume distribution mean the effect is small. 
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All the analysis presented thus far in this sub-section has been based on log-normal 
models. It is useful to consider the errors in /Jcoj produced by uncertainties and noise within 
a single research worker’s data or by uncertainties arising because of contradictions between 
data obtained by different research workers. It is also necessary to consider the fact that 
many measured distributions do not include particles of sufficiently large radii, and such 
distributions must either be discarded, or extrapolated to larger radii. 

Fig. 6.12 shows examples of noise in data obtained by Cress (1980). At the large particle 
end of the size distributions considerable statistical fluctuation is evident. Table 6.2 shows 
values of pcoj calculated for the profiles at heights of 1.6 and 5.6 km, using both the upper 
and lower envelopes of the size distributions. The modeled values for Pco^ differ by more 
than a factor of two in these cases. (Note the values shown in the Table are for extrapolated 
curves; the effect of extrapolation is discussed below.) Another source of uncertainty is 
shown in Fig. 6.13; this shows impactor data obtained, at three different altitudes, for 
stratospheric aerosols by Gras and Ayers (1979) and analyzed using two different tech- 
niques, a third curve in each diagram shows data obtained on the same aerosols using an in 
situ white light technique. All three curves are different at each altitude, particularly c from 
either a or b- Table 6.3 shows the corresponding values for Pcoi, which, at two of the 
altitudes, differ by more than an order of magnitude. 

An even more graphic illustrationn of the uncertainties in /Jcoj is obtained by considering 
the stratospheric data of Hofmann and Rosen which was discussed in Section 3.2.2. In that 
section we saw that Hofmann and Rosen have reported the existence of particles with radii 
between 1 and 2 (xm in greater numbers than had previously been expected. Fig. 6.14 shows 
value of Pco 2 modeled from an example of Hofmann and Rosen’s data. In this example, Pcoj 
has been calculated for: (a) a log-normal distribution fitted to the cumulative number 
distributions obtained at radii of .15 and .25 pm; (b) a slightly truncated log-normal distribu- 
tion with the appropriate number of particles removed and placed in the radius ranges 0.95 
-*• 1.2 pm and 1.2 -*■ 1.8 pm. Corresponding values for Pcn^ are also shown for each section 
of the distribution. It can be seen that inclusion of the large particles increases Pco^ from 2.5 
X 10‘“ m'' sr‘‘ to 55.6 x 10 m'' sr‘‘. Analysis of another example of Hofmann and Rosen’s 

data produced a similar, although not quite as large, increase. 

One problem that occurs is whether it is permissible to use size distribution data that does 
not extend to sufficiently large radii. It is necessary to inquire whether such distributions 
may be extrapolated to larger particle sizes with success. Two methods of extrapolating 
measured size distributions have been considered, which are shown in Fig. 6.15. Since the 
particle size distribution approximates to a double log-normal, it might be expected that a 
log-normal extrapolation would be more accurate than a power law extrapolation. This is 
normally the case and indeed a power law extrapolation can be very unreliable. It has to be 
fitted to the last two or three measured points in the distribution, which as we have seen, can 
be noisy. The resultant extrapolation may lead to very erroneous results. A log-normal ex- 
trapolation is normally more reliable but, if the data are noisy, it may not be possible to fit a 
log-normal. In this case the data are better without extrapolation and possibly should not be 
used at all. Table 6.4 shows the effects of both a log-normal and a linear (power law) ex- 
trapolation on three typical size distributions. The distribution measured by Cress extends to 
a radius of 5 pm. Both a linear and a log-normal extrapolation are possible, neither changes 
the modeled value for Pco^ by a very large factor. The distribution of Gras and Laby is very 
different; it does not extend beyond a radius of 1 f^m. A log-normal extrapolation is not 
possible and the linear extrapolation produces a value of Pco^ that is an order of magnitude 
greater than the unextrapolated value. It is clear that the distribution should not be used, 
either extrapolated or unextrapolated. The third distribution, although extending only to a 
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radius of 1.8 p»m, can be extrapolated by both methods. Again ftco 2 is increased by extrapola- 
tion, but not greatly. In general, it can be said that extrapolation is possible and useful for 
tropospheric aerosol if the measured size distribution extends beyond the peak of the coarse 
particle mode. If the distribution does not extend that far, it should not be used at all. 

6.2.2 /?co 2 calculated from lidar measurements 

Fig. 6.16 shows the results of calculation of the ratios of the backscattering function at 
10.6 /im to the backscattering function at 0.6943 for the same ranges of log-normal 
models and aerosol compositions as used previously. Fig. 6.17 shows similar curves for the 
ratios of the backscattering function at 10.6 (.tm to the backscattering function at 1.06 fim. In 
both cases, the range of variation is somewhat less than that shown in Fig. 6.9 but is, never- 
theless, still about two orders of magnitude over the naturally observed range of log-normal 
modes, even when the curve for soot is excepted. 

A similar range of conversion factors is found when calculations of these ratios are carried 
out for a large number of measured size distributions. The size distributions have been divid- 
ed into three classes — stratospheric aerosol, marine free tropospheric aerosol and continen- 
tal free tropospheric aerosol, and Figs. 6.18 and 6.19 show histograms of the frequency of 
occurrence of a given value for the ratio (calculated for sulfuric acid for the stratospheric 
aerosol, for dust for the other aerosols). This range of variation in the conversion factor has 
been foreshadowed in Section 6.2, where it was observed that the aerosol sizes responsible 
for the bulk of the scattering at 10.6 pm were not the same as those responsible for the bulk 
of the scattering at shorter wavelengths. Even though the data has been sub-divided into the 
classes shown, it is still highly scattered. The stratospheric ratios appear consistent except for 
the fact that these ratios have been calculated assuming no large particles with radii of the 
order of 1 pm. Inclusion of these particles increases the conversion factor by approximately 
one order of magnitude. The marine free tropospheric aerosol, in which the ratio of the 
numbers of large to small aerosols can vary over a very wide range, shows a large variation in 
the conversion factor. Values for the continental free tropospheric aerosols show consisten- 
cy, but this should be used with caution, as many of the values at the peak of the histogram 
come from a single measurement technique (single stage impactor) and may be artificially 
concentrated. 

6.2.3 Pco 2 calculated from extinction and optical depth measurements 

Calculations have been carried out for the ratios of p,o.t to extinction at 1.0 pm in an ex- 
actly analogous manner to those presented in the previous section. The results of these 
calculations are shown in Figs. 6.20 and 6.21. The same comments and conclusions apply to 
these graphs as to the earlier ones and need not be repeated. Fig. 6.22 shows the results of 
calculation of the ratio of the backscattering function at 10.6 pm to the extinction coefficient 
at 10.6 pm. A considerable range of variation (about 1 Vz orders of magnitude) is observed, 
even though both values are calculated at the same wavelengths. 

Detailed calculations of /),o.6 from optical depth measurements have not been attempted. 
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6.2.4 Effects of uncertainties in shape 

Although the smaller sized aerosol particles, especially the liquid phase aerosols, can be 
assumed to be spherical, most naturally occurring scattering particles are not always perfect 
spheres. The shape of aerosol particles can be seen by scanning electron microscope 
photomicrographs or other techniques. An analysis of the samples collected from cascade 
impactors aboard aircraft has shown that the small particles are more nearly spherical in 
shape but are not really smooth spheres and the large particles are irregularly shaped ag- 
glomerates (Wood, 1979). 

In spite of the fact that most natural aerosols are not perfect spheres, a vast amount of 
theoretical calculation has been made for perfect spheres. This is because Mie theory calcula- 
tions are easy to perform and in some cases the qualitative, and even quantitative, 
characteristics of particle scattering are adequately represented by equivalent spheres. For 
most irregularly shaped particles, the calculation of exact backscatter is almost impossible 
even with the use of today’s high speed computers. 

.A recent report by Huffman (1981) shows that the measured volume extinction coeffi- 
cients for AI 2 O 3 are in general slightly larger than those calculated from Mie theory, with the 
exception that, near the absorption line at A = 13.3 /um, the calculated volume extinction 
coefficient is about ten times the measured value. A comparison between the measured 
angular scattring pattern at A = 0.6328 /im for irregular particles with those calculated by 
Mie theory shows that the calculated backscatter is about two to three times the measured 
value (Chylek, 1976). For comparisons at 10.6 /im, the experimentally measured backscatter 
from an ice cloud containing minute crystals (2.3 ^m mean dimension) is only slight smaller 
than that predicted from Mie theory (Sassen, 1981). 

Since the backscattering cross section is a function of the size parameter a = (2nr)/A 
where r and A are radius of the particle and wavelength of radiation, respectively, and the 
size parameter of aerosols at A = 10.6 j^m is equal to only 1/10 of the size parameter of the 
same aerosols at visible wavelengths, we expect the effect of shape on backscatter at 10.6 pim 
should be far less significant than that at visible wavelengths. As mentioned earlier, at the 
visible wavelength A = 0.6328 pim, the difference between the measured backscatter and the 
backscatter predicted by Mie theory is only a factor of about 2. We tentatively conclude 
therefore that the errors produced by Pco 2 due to the irregularity of shape in aerosol particles 
are small in comparison with other sources of error. Only if the wavelength is close to that of 
an absorption line, may the actual backscatter from irregularly shaped particles be 
significantly different from the value predicted by the Mie theory. 
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Table 6,1: Effect of varyinh the width of the log-normal model 


A= 10,5 m . DUST AEROSOL 

(a) fixed number distribution mean, Tq = 0.631 Mtn, 1 PARTICLE PER CH^ 


1.5 

8.73 

X 

.10'^ 

m" 

•1 

sr" 

•1 

1.75 

19.56 

X 

10-9 

m“ 

1 

sr" 

1 

2.00 

31.79 

X 

10'^ 

m“ 

1 

sr“ 

1 

2.25 

46. 9A 

X 

10"^ 

m“ 

1 

sr" 

1 

2.5 

65.14 

X 

10"9 

in“ 

1 

sr" 

1 


3 3 

(b) fixed volume distribution mean> = 2.51 mj 1 per cm 




1.5 

4.10 

X 

10 "^ 


sr' 

-1 

1.75 

3.39 

X 

10 "^ 

m“^ 

sr' 

■1 

2.00 

3.68 

X 

10 "- 


sr" 

•1 

2.25 

3.49 

X 

10 "^ 

m"^ 

sr' 

■1 

2.5 

3.34 

X 

O') 

1 

o 

m"^ 

sr' 

■1 



TABLE 6.2: 


A COMPARISON OF Bcq2 VALUES CALCULATED FROM THE LOWER 

AND UPPER ENVELOPES OF THE REPORTED AEROSOL SIZE DISTRIBUTIONS 

WITH NOISE AT THE LARGE PARTICLE END'^ 


Height 

(tan) 


Bc 02 gr-l)2 

calculated from lower 

envelope calculated from upper envelope 

1.6 

0.816 X 10“® 

1.57 X lO”® 

5.6 

0.671 X lo"® 

1.42 X 10“® 


■^The aerosol size distributions are given by Cress [AFGL-TR-80-0178 (1980) ) 
at 0850 Z on July 7, 1977. (See Fig. 6.12) 

2 

Aerosol size distributions are extrapolated to 30 ym and aerosols are 
assumed to be dusts. 


TABLE 6.3; A COMPARISON OF BcOs VALUES CALCULATED FROM AEROSOL SIZE 
DISTRIBUTIONS DEDUCED WITH DIFFEREOT SIZING TECHNIQUES 


Height 

()cm) 


10 

16 

23 


3c 02 (ra ^ sr ^) 


from Gras & Ayer 
technique 


from Bigg & Ono 
technique 


from In situ white- 
light scatter 


1.568 X lO' 
3.716 X 10 


6.062 X 10 


-11 


■12 


■12 


4.068 X 10 
2.184 X 10 
3.785 X 10 


-12 


-13 


-13 


1.376 X 10 


.-11 


2,866 X 10 


-12 


3.037 X 10 


-12 


1 


Size distribution of aerosol particles are obtained on February 
at Mildura (34.2«>S, 142. 1»E) by Gras and Ayers (1979). Aerosol 
are assumed to he 75% H^SO^, (See pig. 6.13) 


10, 1978 
particles 
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TABLE 6.4: EFFECT OF EXTRAPOLATION ON THE CALCULATED VALUES OF 3 CO 


Data Source 

Maximum Particle 

Altitude 

6c02 

(10-9 n,-l 

sr-1) 



Radius Measured 


No. Extrap. 

Log-normal 

Extrap. 

Linear Extrap. 

Cress, 1980 

5 ym 

6 km 

1.78 

1.93 


2.11 

Gras & Laby 
1981 

1 um 

4-10 km 

0.017 

not 

possible 


0.16 

Patterson et 
1980 

al. 1.8 um 

5 km 

0.60 

0.66 


0.72 



NUMBER 


VOLUflE 



Figure 6. 1 Key parameters of log-normal size distributions and the relationship between 
the number and volume distribution. dN(r)/d log(r) = A„ exp f-k In^ (r/rf]; k = 1/(2 In^ 
oj; r, = r, exp 3(ln if o, = 2.00; r, - 4.23 r,. 



Figure 6.2 Size distributions selected as typical of those occurring in the free troposphere 
(Cress, 1980; Patterson et al., 1980). 
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Figure 6.3 Histograms showing the relative backscatter and extinction, at four different 
wavelengths, from different particle size ranges for the Cress spring average size distribu- 
tion. 
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Figure 6.4 Histograms showing the relative backscatter and extinction, at four different 
wavelengths, from different particle size ranges for the Patterson et al. (1980) continental 
free troposphere size distribution 
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Figure 6.5 Histograms showing the relative backscatter and extinction, at four different 
wavelengths, from different particle size ranges for the Patterson et al. (1980) marine free 
troposphere size distribution. 
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Figure 6. 6 Variation of /?coj with refractive index for the Cress spring average size distribu- 
tion. 
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Figure 6, 7 Variation of ficoi with refractive index for the Patterson et al. (1980) continental 
size distribution. 
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Figure 6.8 Variation of ftco 2 with refractive index for the Patterson et al. (1980) marine free 
troposphere size distribution. 
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Figure 6. 9 Variation, with mode radius of Pco^for constant number density for log-normal 
size distributions and five common aerosol materials. 


</i 

I 


MODE RflDIUSt V) .Mm 

lO-i 10“ 10' 


10' t 


10 “ k 


IC- 


10 ' 


STHTi — r 1 1 min — i i i i!iii| 

? PURE M&7ERUU. 


" _ 

V«NKL£R-NtSSON 

_ 

- 

MDCTURE 

- 

=- 


-<»-= 

I 




^ .'7 o A 

- 


- 

— 


-o“= 

E 0 SULfURC «3D (75%) 


!□ 

wster soluble 

_ 


(SHETTU a FENN) 



AMMONIUM SULPHATE 

I 

- 

1 mml — Ul-iJlIll] — 1_1 

> 1 nil 


10 -“ 10 -' 10 “ 
mode RfiOIUS(N) 


10 ' 

^m 


' Figure 6. 10 Ratios of ftco^ for different aerosol materials to that for dust. This ratio is 
shown as a function of log-normal distribution mode radius. 
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Figure 6.11 [ico^ per unit volume of aerosol for log-normal size distributions, shown as a 
function of number and volume mode radius. 
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Figure 6. 12 Example of aerosol size distribution measurements showing statistical noise at 
the larger radii particles (Cress, 1980). 
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Figure 6.13 Stratospheric aerosol size distributions attained by Gras and Ayers (1979) at 
Mildura, Australia. Distributions marked a and c were determined from impaction samples, 
those marked a using a new sizing model and those marked c the Bigg and Ono (1974) 
model. Distributions marked b were determined independently using in situ white-light scat- 
ter sizing of individual particles. 



Figure 6.14 Contributions to for the Hofmann and Rosen cumulative size distribution 
at an altitude of 13. 7-15.8 km shown in Fig. 3. 15 
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Figure 6. 15 Methods of extrapolating measured size distributions. 
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Figure 6. 16 Backscatter ratios /?,o.6//Jo.69 shown as a function of log-normal mode radius. 
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Figure 6.17 Backseat ter ratios /3io.6//Ji.o« shown as a function of log-normal mode radius. 



Figure 6. 18 Histograms of the frequency of observation of various values of the ratio 

ftio.i/ /Jo-64- 
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Figure 6. 19 Histograms of the frequency of observation of various values of the ratio 
/J| 0 -«/ Pl- 06 - 



Figure 6.20 Back scatter to extinction ratios fito-Jox.m shown as a function of log-normal 
mode radius. 
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Figure 6.21 Histograms of the frequency of observation of various values of the ratio 
/Jl0*6/ 0 \ .00. 



Figure 6.22 Back scatter to extinction ratios /?io.o/o,o.6 shown as a function of log-normal 
mode radius. 
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7. DETERMINATION OF pco. 


7.1 General Behavior 

Fig. 7.1 shows a scattered plot of the calculated values of pcoj for a large number of pub- 
lished size distributions (Bigg, 1975; Grass and Laby, 1981; Farlow et al., 1981; Rosen and 
Hofmann, 1981; Hofmann and Rosen, 1981; Patterson et al., 1980; Cress, 1980; Blifford 
and Ringer, 1969; Blifford, 1970; DeLuisi et al., 1976). Values in the troposphere have been 
calculated assuming a mixed dust/water soluble aerosol and those in the stratosphere have 
been calculated assuming a sulfuric acid/ water aerosol. The figure has been sub-divided into 
regions and it is of particular interest to note that the values for the remote Pacific Ocean fall 
in a Pcoi region quite distinct from that in which the continental free tropospheric values fall. 
(The continental free tropospheric region includes values obtained over the ocean close to 
continental land masses.) It is also important to note that a minimum value of approxi- 
mately 10 “ m ' sr ‘ is found for both the stratosphere and the remote free troposphere. 
Stratospheric values in the Figure are somewhat confused as we have shown the results of 
calculations based on the Hofmann and Rosen data with and without the inclusion of the 
larger particles. Data from Farlow et al. has also been shown with an arrow and two values 
in order to illustrate the effect of the relatively small uncertainty of 15'7o in the mode max- 
imum which is discussed in their paper. Relatively little boundary layer data has been 
included, but the vertical arrows in the bottom right-hand corner of the figure show results 
of calculations based on DeLuisi et al.’s (1976) data, representing vertical averages including 
part of the boundary layer. It may also be noted from the figure that there is a scarcity of 
measured data between 7 km altitude and the tropopause. Also shown on the figure are ex- 
perimental data obtained by MSFC (Bilbro, 1981) and NOAA (Post et al., 1982). This data 
were obtained in June and July 1981. One set of Rosen and Hofmann’s data was obtained in 
July 1981 and it should be noted that the experimental curve passes between the two values 
representing calculations with and without the large particles. Agreement between the ex- 
perimental and the modeled values is good in the upper troposphere and in the boundary 
layer but relatively poor in the lower free troposphere. It is quite likely that much of this 
discrepancy is due to a seasonal variation. Blifford and Ringer (1969) show a very strong 
seasonal variation in the continental aerosol over the USA, particularly in the lower free 
troposphere. This variation has its maximum in June and July when values about five times 
the annual mean were obtained. It is thus likely that the experimental data represents a sum- 
mer maximum; the same is not true for the modeled values. 

Effects due to humidity and temperature changes have not been included in this figure. As 
we shall see in Section 7.5, such effects are normally not of major importance in the free 
troposphere, provided regions of cloud are excluded. 

7.2 Tropospheric Behavior 

In Section 3.2.2 we presented data obtained by Patterson et al. (1980) for the mass aerosol 
loading in the free troposphere over the Pacific Ocean. This was shown in Fig. 3.7 and the 
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right-hand ordinate of the figure shows the corresponding values for [icoi- This has been ob- 
tained by calculating a mean conversion factor (shown in Table 7.1) from mass loading to 
Pco 2 using the size distributions obtained by Patterson et al. (1980) over the Pacific Ocean. 
This factor is undoubtedly variable, probably systematically with latitude and thus the scale 
shown should be taken as approximate only. Nevertheless, it shows clearly the major fact of 
the latitude variation and the probably spread in the minimum value of which is again 
around 10 " m ‘ sr'. 

In Section 3.2.2 we also showed Figs. 3.13 and 3.14 containing the probability distribu- 
tions for two sets of aerosol observations. In Fig. 7.2 these are shown again, converted to the 
equivalent backscattering function. Also superimposed on this figure are the experimental 
probability distributions for pcoi obtained by Post et al. (1982). Points to note are 

(a) The similarity in slopes between the NOAA experimental data and the modeled results 
from Cress 

(b) The lack of height variation in the values of ftcoi based on the Cress data in contrast to 
the NOAA 4 km and 7 km data 

(c) The differences between the values of ftcoj based on the Patterson et al. and the Cress 
continental data reflecting the lack of large particles in the Patterson et al. data 

7.3 Stratospheric Behavior 

Fig. 4.2 shows the long-term stratospheric lidar aerosol variation between 1974 and 1981. 
A mean conversion factor for p,o. 6 /fto .69 has been calculated using stratospheric aerosol data 
without the inclusion of the Hofmann and Rosen large particles (Table 7.1). This has been 
used to generate the pcoj scale shown on the right-hand side of the figure. Once again, a 
minimum value of pcoj of about 10 " m ' sr' is obtained. Shortly after a significant volcanic 
eruption this may increase to 10 '® or 10’ for a period of a few months. (If Hofmann and 
Rosen’s large particles are included in the modeling, even higher values will be obtained.) 
More measurements are required on the variations in the stratospheric particle size distribu- 
tion, particularly for particles with radii ~1 ^m before this figure can be regarded as 
authoritative. 

7.4 Effect of Wavelength Variation 

The backscattering functions have been calculated for CO; wavelengths ranging from A = 
9.1 fim to 11.1 ^m using different aerosol substances and size distributions. The Mie code 
has been used with the refractive indices given in Table 3.5. In the first set of calculations, we 
have assumed the aerosol concentration is one particle per cc and aerosol size distributions 
to be log-normal, with geometric mean standard deviation o, = 2.0 and mode radii varying 
from 0.05 pim to 5.0 /^m. As illustrated in the plots shown in Figs. 7.3, 7.4, and 7.5, for some 
materials, such as sulfuric acid with 75% H2SO4 and 25% H2O, and dust, the variation of 
Pco 2 with CO2 wavelength is small (a factor of less than about 3) regardless of the assumed 
aerosol size distribution. However, for other materials, such as (NH4)2S04, the variation of 
Pco 2 with CO2 wavelength is size dependent. For log-normal aerosol size distributions with 
mode radii at about 0.7 /im, the variation is less than a factor of 2. However, for aerosols 
with large or small radii, it is possible to get an order of magnitude variation. This observa- 
tion can also be seen from the variation of Pco 2 with wavelength, based on measured aerosol 
size distributions (Figs. 7. 6-7. 8). The same three size distributions have been selected for this 
study as were mentioned in Section 6.2. They are, a spring averge continental distribution 
with a large number of the larger particles given by Cress (1980), a continental distribution 
with fewer large particles and a maritime free tropospheric distribution with even fewer 
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larger particles both by Patterson et al. (1980). For the Cress spring average size distrribu- 
tion, the variation of the backscattering function with CO 2 wavelength is always less than a 
factor of 2, regardless of the aerosol substances. For the Patterson et al. size distribution, 
the backscattering functions for water soluble and (NFl 4 ) 2 S 04 aerosols may vary by about an 
order of magnitrude with the CO 2 wavelength, but for other aerosol materials the variation 
is small. 


7.5 Effects of Humidity and Temperature Variations 


As mentioned previously in Section 3.5, a dry but hygroscopic aerosol particle will grow as 
soon as the ambient humidity passes the threshhold value required for deliquescence. When 
aerosols remain in the solution state, their water vapor pressure is in equilibrium with the 
ambient water vapor pressure. The size and effective index of refraction will change follow- 
ing a change in the ambient humidity and/or temperature. Consequently the backscattering 
function and extinction coefficients of the growing or evaporating aerosols change. 

Previously we have developed a numerical model to study th growth and evaporation ef- 
fects on the extinction of 1.0 solar radiation traversing stratospheric sulfuric acid 
aerosols (Yue and Deepak, 1981); we have now extended this model to study both the 
backscatter and extinction, to study not only the wavelength used in the SAGE and SAM II 
satellite experiment (A = 1.0 f^m) but other wavelengths of radiation including A = 10.6 jim 
(CO 2 laser wavelength), A = 1.06 (YAG laser), and A = 0.6943 (Ruby laser). In addition, 
we are able to study not only H 2 SO 4 aerosols, but NaCl and (NH 4 ) 2 S 04 as well. 


We assume that initially all particles have the same composition and their size distribution 
is given by 



l n^(r/r„) -■ 

2(lno*)^ J 


(7.1) 


where A is a constant, and r„ and o, are mode radius and standard deviation, respectively. 


As ambient conditions change, the size of a particle changes from r to r' as discussed in 
Section 3.5. Consequently, the aerosol size distribution changes to n'(r) given by 



l nMr/(fr„)] n 
2(lno,)^ J 


where f is the growth factor defined in Section 3.5. 


(7.2) 


The percentage change of backscattering function when ambient conditions change is then 
given by 


P = ^ X 100% 

P 


a t 


ln^[r/(frJl 

2(lno,)^ 


}dr- T 

.r2 

J r. A exp J - I dr 

r 2nno.')^ J 


dr 

2(lno,)^ J 


X 100% 


Some of the percentage backscatter changes, plotted as a function of temperature or 
relative humidity change, are presented in Figs. 7.9 to 7.12. These figures illustrate that the 
change of backscattering function due to a change in ambient conditions is negligibly small 
at low relative humidities. However, when the relative humidity is 90% or higher, the change 
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in the backscattering function may reach more than 100*7o for a few degrees change in 
temperature or a few percent change in relative humidity. The plotted curves also 
demonstrated that the change in backscattering function is size dependent; the smaller the 
aerosol particles, the higher the backscatter change corresponding to the same ambient con- 
dition change. 

In order to see whether high relative humidity actually occurs in the atmosphere, we have 
studied the relative humidities reported by Cress (1980). The result is plotted in Fig. 7.13 
where each dot in the height versus relative humidity plane represents a value reported by 
Cress. Also plotted in Fig. 7.13 are the median for Cress’ data and the profiles of humidity 
in January and July reported by Newell et al. (1972). It can be seen that the higher the 
altitude, the less the chance of occurrence of high humidity. At high altitudes there still ex- 
ists some instances where the relative humidity is close to 100®7o. Consequently, the effect of 
humidity on backscatter should not be ignored, especially for aerosols near the ground. 

The percentage ocurrence of humidity at 1 km is plotted in the form of histogram in Fig. 
7.14. Also plotted in Fig. 7.14 are the ratios of backscattering function at a given relative 
humidity to the backscattering function at 0% relative humidity, using the Cress spring 
average aerosol size distribution with two aerosol materials: (NH4)2S04 and NaCl. The deli- 
quescence starts at a relative humidity equal to about 77*7o and 82% for NaCl and 
(NH4)2S04, respectively. The ratio of the backscattering functions for (NH4)2S04 drops 
down to less than 1 before it increases and gradually reaches values greater than 1 at a 
relative humidity greater than 95%. This surprising phenomenon is due to the fact that the 
refractive index of dry (NH4)2S04 is about 1 .98 but that of H2O is only about 1.18. Although 
the absorption of water after passing the deliquescent point increases the cross section of the 
particles, the decrease in refractive index plays a more dominant role, resulting in a decrease 
in backscattering function. The behavior of NaCl aerosols is less unexpected, showing a 
large change in backscattering function following the same change in relative humidity in a 
humid environment. , 
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TABLE 7 . 1 ; CONVERSION FACTORS USED TO OBTAIN < 3 ^ 0 ? 


Atmospheric Region 

MARINE FREE TROPOPAUSE 
AEROSOL 


Conversion Factor 

FOR 1 #xGM'^ 

= 2.0 X IO'^^m'^SR’^ 


STRATOSPHERIC AEROSOL ^10.6/^0.69 

= 3.0 X lO'^ 

^ 10 . 6 ^^ 1.00 
= 1.0 X io’^^sr"^ 

N O TE 

AEROSOL DENSITY IN /iG.M'^ = AEROSOL VOLUME CONCENTRATION 

IN mM^.CM’^X P 
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BfiCKSCFiTTERING FUNCTION, 


Figure 7. 1 Modeled backseat tering functions as a function of altitude for a large number 
of particle size distributions (see text for full references). Also shown are two curves for ex- 
perimental determined values of ftcoi (Bilbro, 1982; Post et al.. 1982). 
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Figure 7.2 Cumulative probability distribution for modeled (ico^ derived from the data 
shown in Figs. 3. 13 and 3. 14. Also shown are the cumulative probability distributions ob- 
tained by NO A A (Post et al., 1982). 
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MODE RfiDIUS. 

Figure 7.3 Pco, versus mode radius for H2SO4 aerosols and different CO2 wavelengths. 
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Figure 7.4 pco2 versus mode radius for dust aerosols and different CO2 wavelengths. 
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Figure 7.5 ftcoi versus mode radius for (NHt)iSOt aerosols and different COi wavelengths. 



WRVE LENGTH, 


Figure 7.6 Pro, versus wavelength for the spring average aerosol size distribution reported 
by Cress (1980) and for different aerosol compositions. 
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NAVE LENGTH, 

Figure 7. 7 fico. versus wavelength for continental aerosol size distributions reported by 
Patterson et al. (1980) and for different compositions. 



Figure 7.8 Pco^ versus wavelength for marine aerosol size distributions reported by Patter- 
son et al. (1980) and for different aerosol compositions. 
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Figure 7.9 Percentage ftco 2 change versus temperature change for HiSOt aerosols, using 
log-normal size distributions with mode radii equal to 0. 05, 0. 5, and 5. 0 yon, at temperatures 
15°C and 25°C and at a relative humidity of 20%. 



Figure 7 . 10 Percentage /Jcoj change versus relative humidity change for HiSO* aerosols, 
using log-normal size distributions with mode radii equal to 0.05, 0.5, and 5.0 yon, at 
temperatures 15°C and 25°C, and at a relative humidity of 60%. 
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Figure 7.11 Percentage Pcot change versus percentage relative humidity change for NaCl 
aerosols, using log-normal size distributions with mode radii equal to 0.05, 0.5, and 5.0 jjtm, 
at temperatures 15°C and 25°C, and at a relative humidity of 90%. 



Figure 7. 12 Percentage (icoi change versus temperature change for (NHfiSOt aerosols 
using log normal size distributions with mode radii equal to 0.05, 0.5, and 5.0 yon, at 
temperatures 15°C and 25°C, and a relative humidity of 90Wo. 
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Figure 7. 13 The humidities at different heights as reported by Cress (1980) and Newell 
(1972). Each dot represents an ambient condition reported by Cress and the solid line is the 
median value for this data. The short- and long-dashed lines are the 40°N January and July 
values presented by Newell et al. (1972), respectively. 
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Figure 7.14 Histogram of the percentage of cases of a given relative humidity at an altitude 
of 1 km, reported by Cress (1980). The curves represent the ratio of ftco. at a given relative 
humidity to ftcoj at 0% relative humidity. O is for (NHJ^SO^ aerosols and □ is for NaCl 
aerosols. 
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8. CONCLUSIONS 


Our conclusions from this study have been assembled in the form of tables, each covering 
a separate aspect of the work. The division is as follows. 

TABLE 8.1: Deduction of Pcoifrom Aerosol Measurements. This table summarizes the 
main facts to emerge concerning the problems of modeling /?co,, using measurements of 
aerosol particle size and refractive index as a basis. 

TABLE 8.2: Main Features of Lidar and Extinction Measurements of Aerosols in the 
Lower Atmosphere. This table summarizes the main characteristics of the available data set 
of lidar and extinction measurements made on atmospheric aerosols. 

TABLE 8.3: Deduction of (ico^from Backscattering Measurements at Other Wavelengths 
and from Extinction Measurements. This table summarizes the problems of determining /Jcoj 
using other optical measurements as a basis. 

TABLE 8.4: General Remarks Concerning the Modeling of Pco^- This table summarizes 
the range of accuracy achievable when modeling is done. It lists the major information gaps 
and makes recommendations for filling these gaps. 

TABLE 8.5: Effects of COi Wavelength and Ambient Condition Changes on pco 2 - This 
table summarizes the changes in the backscattering function due to a change in radiation 
wavelength from 9.1 /^m to 11.1 jim. and identifies parameters affecting the change in pcoi 
when the ambient conditions change. 

TABLE 8.6: Sources of Error in the Calculation of pcoi- This table lists (separately from 
what may already have been shown in the other tables) the possible sources of error in the 
calculation of pcoi and the magnitude of their effect on /Icoj- 

TABLE 8.7: Modeled Values of Pco 2 In the Atmosphere. This table summarizes the infor- 
mation that we have been able to obtain concerning the actual behavior of Pco 2 in the at- 
mosphere, its magnitude, its variation, and the major information gaps and problems. 
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TABLE 8.1 DEDUCTION OP 6^02 FROM AEROSOL MEASUREMENTS 

1. At a wavelength of 10.6 ym, particles with radii greater than 0.5 ym 
are responsible for 50-95% of the scattering from tropospheric aerosols. 

2 . 6 CO 2 For 1 particle per cm varies from 10 ^ 10 ^ m ^ sr ^ for 

naturally occurring aerosols, the lowest values being associated with 
the smaller particles. 

3 . for an aerosol volume concentration of 1 ym^ per cm varies 
^ —10 —8 “1 “1 

from 10 to 10 m sr for naturally occurring aerosols. In 
certain circumstances, analysis in terms of aerosol mass distribution 
rather than mmiber distribution may be preferred. 

4. Shape is not an important factor in the calculation of 3 c 02 ' except 
possibly near an absorption line. (It may affect extrapolation 
from shorter wavelengths.) 

5. Accurate calculation of 3 c 02 J^equires accurate knowledge of dn/d log r 
to at least 1 ym and probably 2 ym for stratospheric aerosols and to 
at least 5 ym and probably 10 ym for tropospheric aerosols. 

6 . Composition is an important factor only for aerosols with a number 

mode radius greater than 1 ym. For a mode radius of 5 ym the difference 
in 6 (-q 2 between a water soluble and a non-water soluble aerosol may 
approach an order of magnitude. 

7. Extrapolation of aerosol size distributions to large sizes may lead 
to significant errors. 
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TABLE 8.2 LIDAR AND EXTINCTION MEASUREMENTS 


I. Lidar: 

1. Mainly confined to U.S.A., Europe and Japan 

2. Principal wavelengths are 0.6943 ym, 1.06 ym (and 10.6 ym) 

3. Long term data sets exist for the stratospheric aerosol 

4. Numerous boundary layer measurements exist 

5. Due to normalization problems, data for the free troposphere 
(and particularly its upper levels) is of limited value 

6. Recent measurements of aerosols, using airborne lidar systems, 
at both visible and infrared wavelengths, show that this is a 
powerful technique for investigating many aerosol characteristics. 

II. Extinction: 

1. The most significant extinction data set, with global coverage, 
is that obtained at a wavelength of 1 ym from the SAGE and SAM II 
satellite systems. 

2. Although the majority of data obtained by SAGE and SAM II is 
confined to stratospheric regions below 30 km, under favorable 
conditions, extinction data to the earth's surface can also be 
obtained. 
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TABLE 8.3 DEDUCTION OF 6c02 BACK- SCATTERING MEASUREMENTS AT 

OTHER WAVELENGTHS AND FROM EXTINCTION MEASUREMENTS 

1. Particle sizes responsible for most extinction and scattering 

at shorter wavelengths are not the same as those responsible for 
scattering at 10.6 ym. 

2. Values for the ratios of 8 c 02 ^0.69' 06’^^'^ *^1.00 vary 

over 1.5 -»■ 2 orders of magnitude. Some sub-division into aerosol 
class may be made, lowest values of the ratio occur for stratospheric 
aerosols . 

3. Stratospheric lidar measurements may be used to give an approximate 

time history for 6 qo 2 stratosphere. Minimum values may be 

fixed fairly accurately. 

4. The ratio 8^02 ^C 02 varies over 1.5 -*■ 2 orders of magnitude, 

lowest values occurring for stratospheric aerosols. 
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TABLE 8.4 GENERAL REMARKS CONCERNING THE MODELING OF &CO 2 

1. The qualitative behavior of 3^02 modeled including, 

in particular the height variation and major geographical 
variations . 

2. Given present published information on aerosols, errors are half 
to one order of magnitude, even for the best information. 

3. Major information gaps are 

1. Composition and size distribution of large particles 
at all levels above the boundary layer 

2. Seasonal and geographical variations of the aerosol 

3. Aerosols properties at altitudes of 7-12 km 

4. Humidity effects on large aerosols 

4. Considerable variation in the parameters as measured by different 
instruments now exists 

5. Extrapolation of optical measurements at other wavelengths to 1.06 ym 
can lead to errors at least as large as an order of magnitude. 

The error is dependent upon aerosol size distribution and may be reduced 
if this is known more exactly. 

6. It appears desirable to 

1. Make simultaneous measurements for comparison purposes, using 
different sensors in the free troposphere 

2. Attempt to fill the information gaps noted above 

3. Make simultaneous measurements of CO 2 back-scatter and aerosol 
properties to confirm and improve modeling. 
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TABLE 8.5 


EFFECTS OF CO^ WAVELENGTHS AND AMBIENT CONDITION CHANGES ON BACKSCATTERING FUNCTION 


The variation of 6 , the backscattering function at CO wavelengths (9.1 to 

CO 2 2 

11.0 um) depends on the aerosol composition. 

1. For aerosol particles other than water soluble and ammonium sulfate, 

the variation in is less than a factor of 2. 

CO 

2 

2. For water soluble and ammoni\im sulfate aerosols, there may be a one 

order of magnitude change in . 


Effect's of changes in relative humidity and temperature on depends on 

several parameters: 

• Aerosol composition — changes in B from H.,SO., (NHJ_S0^ and NaCl 

CO^ 2 4 424 


aerosols are different. 

Aerosol mode radius — the smaller the mode radius, the larger the change 
“ ®C0 ■ 


Ambient temperature — the lower the ambient temperature, the larger the 

change in B^q . 

^2 

Ambient water content — the higher the ambient water content, the larger 
the change in B„ . However, aerosol deliquescence must occur before 


any change can be observed. 

Change in ambient condition — increase in ambient water content or 
decrease in temperature will in general increase B^^q / except near the 
point of deliquescence. 
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TABLE 8.6 


SOURCES OF ERROR IN THE CALCULATION OF 6c02 


Uncertainties in 

1. Aerosol size distribution and concentration 

— Lack of detailed knowledge of particles with radii > 1 ym make 3^^ 
calculations totally invalid. Even using the best information 
available on these particles, order of magnitude errors may arise. 

2. Aerosol refractive index 

— Except for materials with very high or low refractive indices 
or near absorption lines, errors are not more than a factor of 
two or three. 

3. Aerosol shape 

— Not an important factor (errors ^ a factor of two) except possibly 
near an absorption line. 

4. Measurement uncertainties 

— Sampling problems arise here, as well as uncertainties in the 
technique itself. Order of magnitude errors may occur. 

5. Natural variations 

-- The altitude variation of aerosol concentration and composition 
is reasonably we 11- documented from 0-6 km and in the stratosphere. 
Geographical, seasonal and meteorological variations are not well- 
documented but concentration variations of aerosol orders of 
magnitude definitely occur. 
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TABLE 8.7 MODELED VALUES OF 6 , 


IN THE ATMOSPHERE 


CO2 


1 . 


3 co^ decreases from 10 m sr at the earth's surface to 10 
m sr in the stratosphere. 


2. Values determined from different data sets may differ almost as much 
as the range within a single data set, 

3. There is a scarcity of data for heights between 7 km and the tropopause. 

4. 6 c 02 stratospheric aerosols is very dependent upon the presence 

of large particles with radii > 1.0 ym. 

5. Significant differences exist between values for 6 CO 2 free 

troposphere over the land and over the remote ocean, the latter 
being an order of magnitude less, on the average. 


6 , 


The 10% and 90% probability levels for 
lie 1 1 - 1/2 decades apart. 



in the free troposphere 


7, The height variation of 6 co 2 probability distribution 

of 6 co 2 ^ 9 ^®® reasonably well with measured values. 
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